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ABSTRACT 

This work constitutes a comprehensive investigation of the application 
of cryogenic technology to coils suitable for MHD applications. Detailed 
background studies are conducted in the specific areas of resistive behavior 
of metals, coil structural requirements, cryogenic refrigerator character- 
istics, properties of coolants, and cooling system design. These studies 
were utilized to develope a series of characteristic curves, suitable for use 
in the preliminary design phase of this type of cryogenic magnet. The range 
of the investigation covers coils suitable for generating fields of up to 100 KG, 
for use with MHD channels of cross section between . 01 and 10 square meters. 

Consideration is given to the optimization of operating temperature and 
coil geometry. A simplified procedure for preliminary coil analysis, uti- 
lizing the characteristic curves, is presented in the form of a flow diagram. 
Example problems indicate that "minimum loss" coils appear highly attrac- 
tive for use with seeded-gas MHD generators of 50 Mw output, or larger. 

In general, "optimized" cryogenic coils should provide at least an order of 
magnitude reduction in total power requirements, compared to "conventional" 
copper coils. A reasonable evaluation of present data indicates that event- 
ually, this may be improved by yet another factor of four. 
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I. INTRODUCTION 



A. Background : 

Recent developments in the field of magnetohydrodynamics point to- 
ward an increasing need for high-field electromagnets combining large 
physical size with low power requirements: In particular, the application 
of seeded-gas MHD generators to commercial power generation (1,2) in- 
dicates a probable need for fields of the order of 20-50 kilogauss, with 
channel lengths ranging between 1 and 30 meters and widths between 0. 1 
and 5 meters (the larger sizes corresponding to the lower fields). Further, 
recent interest in the application of MHD to marine propulsion (3, 4) sug- 
gests a possible need for fields of up to 200 KG, with channel dimensions 
scaling upwards from the larger figures cited above. It has been shown 
(2, 5) that the unprecedented size and high field requirements, coupled with 
the demand for operating economy, preclude the use of conventional 
electromagnets for these applications. 

The startling advance of superconductor technology within the last 
decade has led to the consideration of superconducting magnets, to the 
exclusion of all others, for applications such as the above. Indeed, in- 
vestigation of the Vanadium-Gallium system indicates the possibility that 
"hard" superconductors may be capable of operation in fields up to 500 
KG (6). Despite this promise, however, superconducting magnets exhibit 
several less -than-desirable characteristics which have hampered their 
development, to date: 

1. The elements utilized in superconducting alloys and compounds 
are relatively rare, and correspondingly expensive. 

2. The mechanical properties of these alloys and compounds most 

notably, brittleness lead to severe problems in the manufacture and 

9 







1 













operation of high-field coils. 

3. The present "standard of the industry" ten mil wire poses 

serious problems regarding the assembly of physically large coils. 

4. Due to the fact that hard superconductors become relatively poor 
conductors when quenched, the magnetic field of a suddenly quenched super- 
conducting magnet is prone to rapid collapse. This leads, at best, to the 
problem of rapid disposal of a great deal of stored energy, ana can lead, 

at worst, to local arcing and/or destruction of the coil. Subsequently, 

local impairment of the cooling system associated with such a coil if 

sufficient to permit a local temperature fluctuation to the vicinity of the 

quench point could necessitate immediate and total shutdown. This, 

in itself, could be disastrous, particularly in marine propulsion appli- 
cations. 

5. To date, the current density attainable in large superconducting 
coils has fallen far below that predicted on the basis of short-sample 
testing, thus adding an additional complication to the design of high-field 
magnets. 

Although significant advances have been made in the search for solu- 
tions to these problems particularly, as regards coil protection (7), 

and determination of the factors underlying degraded performance in large 
coils (8-14) the successful synthesis of a large, high-field supercon- 

ducting coil has proven to be an elusive, and as yet unrealized, goal. This 
fact has become a critical roadblock to the advance of MHD technology: 
so much so, that background studies leading to the development of a com- 
mercial MHD generator have perforce been conducted with a large, pulse - 
operated conventional copper-coil electromagnet (15, 16). 

An alternative to the utilization of superconducting magnets for the 
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above applications would be the use of so-called "cryogenic" magnets: 
that is, magnets utilizing coils wound of any of several highly-purified 
and carefully annealed metals, designed to take advantage of the marked 
resistance decrease of "ordinary" conductors at cryogenic temperatures, 
as contrasted with superconductivity. Such magnets share the promise of 
relatively economical generation of intense magnetic fields over large 
volumes, but not the problems cited above. Specifically: 

1. Initial material cost is less. 

2. Mechanical properties tend toward ductility. 

3. No significant restrictions need be placed on conductor size. 

4. The metals contemplated are those considered "good" conductors, 
even at room temperature. Thus, resistivity is not critically sensitive 
to minor temperature fluctuations. 

5. No serious scaling problems are known to exist. In fact, many of 
the problems associated with cryogenic magnet design become less severe 
as size is increased. 

Of course, these advantages are offset by relative disadvantages: For 
example, since cryogenic coils do have non-zero resistivity, ohmic 
heating occures within the coil dewar, and an appreciably larger refriger- 
ation system may be required for a given cryogenic coil than for a com- 
parable superconducting coil. This implies both higher operating costs 
and higher equipment costs, with the relative increase depending primarily 
on coil size. Again, larger coils would appear more competitive on a 
relative basis. However, two facts seem pertinent, here: First, even 
relatively small cryogenic coils can be designed to operate considerably 
more economically than comparable room -temperature magnets. Second, 
cryogenic magnets have been built and successfully tested in sizes larger 

than those yet attainable by superconducting technology. 
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B. State-of-the-art: 



Interestingly enough, although at least one study indicates that large 
cryogenic coils might closely approach the anticipated efficiency of the 
yet elusive superconducting coils for large-scale MHD applications (2), 
none have been built for this purpose, and in fact, no detailed study of 
such an application has been published to date. Perhaps this can be 
attributed to the glamour of superconducting magnets. However, con- 
siderable interest has been exhibited in the development of cryogenic 
magnets in geometries suitable for controlled fusion experiments, as 
can be seen from a brief description of coils already constructed in this 
country: 

Los Alamos Scientific Laboratory (A. E. C. )(17, 18 ): 

Size: 2-1/2" I. D. x 7-1/2" O. D. x 5" long. Construction: "Jelly roll": 
120 turns of 5" x . 010" copper foil, separated by inert spacers for 
coolant channels. Coolant: Liquid hydrogen (immersion). Field: 60 KG 
steady, 62 KG for 5 second pulse. Status: Operational. 

Project Sherwood (N. B. S. ) (19,20) : 

Size: 3" I. D. x 12" O. D. x 7" long. Construction: "Tape wound": 14 
stacked pancakes, each 0. 4" thick, wound of 14, 000 turns of 4 mil alu- 
minum foil (insulated by 0. 17 mil paper strip) and encapsulated in epoxy. 
Radial slots machined in inert spacers between pancakes for coolant flow. 
Coolant: Liquid hydrogen (forced flow). Field: Tested to 95 KG. Status : 
Structural failure at 95 KG: magnet burst when Lorentz forces exceeded 
inherent strength of epoxy-bonded A1 foil. 

Lewis Research Center (N. A. S. A. ) (21, 22) : 

Size: 30 cm. I. D. x 90 cm. O. D. x 60 cm. long. Construction: "Jelly 
roll": 12 stacked pancakes, each with 75 turns of 5 cm. x 0. 2 cm. alum- 
inum strip, separated for coolant flow. Coolant: Liquid neon (immersion). 
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Field: Designed for 200 KG. Status : Undergoing preliminary testing. 

Lawrence Radiation Laboratory (U. C. ) (23, 24) : 

Size: 7-1/8" I. D. x 16" O. D. x 12" long. Construction : 77 turns of 
0. 625" square conductor (sodium cast into . 010" stainless tubing), wound 
in four layers with coolant grooves and stainless support cylinders between 
adjacent layers. Coolant : Gaseous helium (forced flow). Field: Designed 
for 100 KG. Status : Failed during testing at 21 KG when Lorentz force 
extruded sodium through small crack in stainless jacket. Being rebuilt. 

The above coils share many common features: All are cylindrical 
solenoids, such as might be suitable for construction of a magnetic mirror. 
All have what has hitherto been considered a large bore, for high-field 
magnets. All but the last (designed by Taylor and Post) are of extremely 
simple design, and even the last design appears simple contrasted with 
what might be required to perform the same task with "conventional" 
electromagnets. 

Taylor and Post are evolving what might be considered the second 

generation of cryogenic magnets: Their design represents the first 

attempt to obtain the best possible performance (25-27) from a cryogenic 

magnet, as contrasted with previous attempts to obtain performance, 

period. They have attempted to optimize their design in two respects: 

first, they have attempted to identify that operating temperature which 

2 

minimizes total power expended in removing I R losses from the coil, 
and second, they have identified the particular geometry which maximizes 
magnetic energy density per unit dissipation for a cylindrical solenoid. 

C. Application to MHD : 

The coil model which will be used throughout this analysis is the 

overlapping cylinder or "saddle coil" model discussed by Stekly, et al. 

13 



(2) (Fig. 1). Conceptually, this model consists of two identical, parallel, 
circular cylinders, carrying a uniform, equal axial current in opposite 
directions, which have been brought together until they actually overlap. 
Since, by superposition, no current flows in the region of overlap, the 
material in this region can be removed to provide the "working area" 
of the coil. Such a configuration confers many advantages for trans- 
verse-field MHD applications: First, referring to the axes depicted in 

Fig. 1, it can be shown that in the working area, the magnetic field pro- 
duced by the overlapping cylinders is uniform everywhere and is given by: 

B x ’ B y =0 ' B z'/Ve a (1) 

where is the effective current density in either cylinder and a is the 

offset between the center of either cylinder and the center of the working 
area. Next, insofar as MHD channels tend to be relatively long and 
narrow, the saddle coil configuration offers a relatively efficient utili- 
zation of conductor material, in that the coil windings are so arranged 
that the individual conductors are kept close to the MHD channel which 
they serve. Finally, the configuration is practical from an engineering 
viewpoint: The indicated cylinders would, in practice, be comprised of 
a large number of individual conductors, so that the "crescent moon" shape 
of the coil lobes could be readily approximated by the assembly of modular 
units. In short, it seems reasonable to assume that actual MHD coils 
would closely resemble such a model, with the degree of similarity prob- 
ably controlled by a compromise between desired field homogeneity and 
ease of assembly. 

It is apparent that many of the characteristics of a particular coil 
will be determined by the application for which it is intended. For example, 
it has been shown (28, 29) that the length of a transverse-field MHD channel 
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FIGURE 1 



Saddle Coil: Side View 
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should be 5-10 times the mean channel width, to minimize degrading end 
effects. It will be assumed here that the contemplated MHD channel is of 
rectangular cross section, with width 2w and height 2h. The largest such 
rectangle which can be inserted into a given working area has: 

-3a + Va 2 + 8R 2 

w = (2) 

4 

and 

A , = 4wV R 2 - (a + w) 2 (3) 

ch 

where R is the radius of either cylinder. The length of the coil associated 
with such a channel will be taken to be: 

(4) 

The total cross section (both lobes) of the coil is related to the above 
parameters by: 

A c] _ = 2R 2 [7T- (2 /- sin 95 ]. 7 - cos *(a/R) (5) 

One additional concept must be included in the discussion of such 
coils: At either end of the coil, the individual conductors must pass up 
and over (or down and under) the MHD channel. It will be assumed here 
that these "end turns" do not contribute to the useful magnetic field, but 
do contribute to total loss. Obviously, one would desire to minimize the 
mean length of the end turns, within the constraints imposed by the above 
parameters. Considering this, the mean length of an end turn will be 
estimated as one half the circumference of the circumscribing circle: 
that is, 

L e ^7T(R + a). 

It is then useful to define a form factor, F, as follows: 

F = 1 + L /L (6) 

e 

Thus, geometrically similar coils can be described by the common value 

of their form factor. Further, total coil volume, etc. , can be estimated 
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as: 



V. = F A ,L 
tot cl 



(7) 



The purpose of this work, then, is to examine the application of 
cryogenic technology to coils such as those described above. The basic 
philosophy of this investigation is to seek that information required to 
support a realistic appraisal of the desirability of cryogenic magnets for 
a particular MHD application, to identify those parameters which can be 
used most conveniently to define a coil during the preliminary design 
process, and to illuminate, in detail, the relationships between those 
parameters and the resultant coil. The range of the investigation covers 
coils suitable for generating fields of up to 100 KG, for use with MHD 
channels of cross section between . 01 and 10 square meters. It might 
be noted that, aside from the apparent geometric differences, such coils 
differ from the solenoidal coils discussed earlier in that the constraints 
imposed imply approximate coincidence in size between the smallest 
coils considered here, and the largest cryogenic solenoids built to date. 
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II. RESISTIVITY OF METALS 



A. Background : 

It is appropriate to begin this analysis of cryogenic magnets with a 
brief resume of the physical reasons which favor the application of cryo- 
genic technology to magnet design: Essentially, the resistivity of pure 
metals is a sufficiently strong function of temperature that refrigeration 
to the cryogenic range can reduce magnet power consumption by more 
than the additional power required to run the refrigerator. For this 
reason, cryogenic magnets can be made to operate with significantly less 
power consumption than comparable "conventional" magnets. 

MacDonald (30) provides an excellent summary of the development 
of the underlying theory: Briefly, Mathiessen (31) first recognized that 
the resistance of metals could be separated into two components, one 
dependent upon the degree of purity, and the other upon the temperature. 

He postulated that these components were independent, and therefore 
directly additive. They have come to be generally known as "residual" 
or "impurity" resistance and "intrinsic" or "ideal" resistance, respect- 
ively. The former is attributed to scattering of conduction electrons by 
local deformation of the crystal lattice (e. g. , by the presence of impurity 
atoms, or dislocations), and the latter, to scattering by thermal vibration 
of the lattice. Thus, the "ideal" resistance goes to zero as the absolute 
temperature goes to zero. 

No significant advance was made in quantitative analysis of the thermal 
behavior of ideal resistance until quantum theory became fairly well 
developed: In 1929, Bloch (32), utilizing a lattice model first postulated 
by Debye (33) in 1911, evolved the first strictly theoretical treatment to 
predict accurately the low-temperature resistive behavior of "ideal" 
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nretals. Bloch's Law states: 

*' T 5 



P = 




x" dx _ h Vm 

(e x -1) (1 - e x ) k 



where ©, the "Debye temperature", is a constant, characteristic of a 
particular lattice. (Actually, Bloch's original solution considered only 
the low- and high-temperature limiting cases, but subsequent investigators 
have found that the above integral expression, derived by him, is essen- 
tially valid for all temperatures. For T «©, the above expression reduces 
top - (T/©) 5 , while for T^©, p~T /©. ) 

Within a few years, Gruneisen (34) noted that Bloch's theory could 
be modified to provide a more accurate prediction of the resistive behavior 
of real metals by use of an empirical correction factor: He suggested that 
the resistivity derived from eqn. (8) be multiplied by a factor of the form 
(1 +Q£T +CX2T ), where and were to be determined experimentally 
for the various metals. While admitting that this technique could not be 
defended on strictly theoretical grounds, he advanced heuristic arguments 
defending the form of the indicated correction, and observed that in most 
cases, the resultant "fit" to experimental data was accurate to within 
about 10% of the value of the residual resistance. (Later investigators 
(35) have introduced an equivalent, though somewhat less convenient, 
concept of describing "real" resistive behavior by use of a variable 
Debye temperature. ) 

Concurrently with the above developments, several investigators 
became interested in the fact that the resistance of metallic conductors 
increases in the presence of a magnetic field. In 1938, Kohler (36) was 
able to show that: 



T" = H 7 



H_ 

r 



(9) 



where r Q = r(0, T) is the total resistance at temperature T, in the absence 
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of any magnetic field, and Ar = r(H, T) - r(0, T) is the resistance increase 
in the presence of the field. Kohler was unable to obtain a general sol- 
ution for the functional form of though he predicted that in the limit 
of small arguments, 

¥(H/r )— a(H/r ) 2 , 
o o 

while in the limit of large arguments, 

^(H/r o )-^b(H/r o )°, 

i. e. , the magnetoresistance (Ar) saturates. Obviously, magnetoresistance 
is more important at cryogenic temperatures than at "ordinary" temper- 
atures, since the reduced value of r Q leads to a more rapid increase of 
magnetoresistance with field for any given metal. Further, though Kohler 
did not obtain a general solution for\j/, mere recognition that the functional 
relationship exists suggests a plotting technique by which the results of 
limited testing may be generalized to predict the magnetoresistance of 
samples of different purity, or at different temperatures. Plots of 
magnetoresistance vs. H/r Q have come to be known as "Kohler diagrams". 
(Equivalently, one can plot magnetoresistance vs. the ratio of mean free 
path in the lattice to the radius of the orbit of a free electron in the mag- 
netic field: However, though this latter concept may be more elegant in 
a physical sense, the former technique is generally more convenient for 
the experimenter. ) 

In 1940, Justi(37) recommended the use of a so-called reduced Kohler 

diagram, in which magnetoresistance is plotted vs. H(r~/ r ), rather than 

y o 

H/r Q . It was his hope that normalizing the field-free resistivity to its 

value at the Debye temperature for the particular lattice would eliminate 

differences due to the variation in "intrinsic" scattering between different 

lattice structures. Unfortunately, such normalization does not cause the 

data, for all metals to fall on a single curve, though the curves for the 
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various-metals, when so plotted, do tend to fall in a rather broad band. 
Nevertheless, many experimenters have adopted the practice of plotting 
magnetoresistance data with field-free resistivity normalized in this 
manner. Actually, little is gained, and much is lost, through this tech- 
nique: In general, authorities cannot agree on the exact value of the 
Debye temperature for a particular metal. Further, since Debye temp- 
eratures for most metals fall well out of the range in which detailed 
experimental resistivity figures are available, the normalizing standard 
for a particular metal must often be generated from the theory. Aside 
from the approzimations inherent in the theory, there is the problem of 
which version to use: Bloch's original version, the Gruneisen modification 
thereto, or some more recent derivation. Thus, one is faced with data 
normalized to a rather fuzzily defined standard, and all too often, pub- 
lished data lack any indication of the normalizing standard used. Even 
when the standard can be determined, this practice obscures any direct 
comparison between different metals. About all that can be said to 
recommend the practice is that it reduces the size of the numbers involved. 

Despite continued interest in the theory of the magnetoresistive 
behavior of metals, little significant advance has been made in the last 
several years. To quote from a second article by MacDonald (38): 

" it appears that the theory is quite unable to give any 

detailed explanation of such bizarre behavior: Wilson says 'It 
does not seem likely that any model simple enough to be tractable 
theoretically would give a magnetoresistance curve of the com- 
plexity of those actually observed. '" 

B. Limitations of the theory: 

The above relations, taken together, comprise the theoretical moti- 
vation for the use of cryogenic magnets: In theory, careful purification 
and annealing will result in a residual resistance several orders smaller 

than the~room temperature resistance of a given metal. Cooling to cryo- 
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geni-c temperatures can make the ideal resistance arbitrarily small 
compared to the residual resistance. Then, if the magnetoresistance 

saturates at a sufficiently low value or, at least, does not rise too 

rapidly with increasing field cryogenic magnets can be made to operate 

much more economically than conventional electroniagnets. 1 

Unfortunately, as might be inferred from the above discussion, the 
approximations inherent in the foregoing concepts preclude their direct 
application to magnet design: 

1. The theory does not quantize the relationship between degree of 
impurity, or extent of cold work, and the magnitude of the resultant 
residual resistance. Nor is anything said regarding the feasibility of 
attaining a specified degree of purity in a particular metal. 

2. Experimental data indicate that Mathiessen's rule is an best an 
approximation. Assuming its validity, one can compute the "ideal" 
resistance-temperature curve for a given metal by extrapolating an 
experimental curve to 0° K to determine the residual resistance, then 
subtracting that constant from the original curve. In general, curves so 
computed from observation of samples of lesser purity (34, 39) lie well 
above those computed from measurements on highly pure samples (40-43), 
at least for the liquid hydrogen-liquid helium range. Maimoni's date 

(48) indicate a similar discrepancy for aluminum samples of the same 
purity, in which the resistance ratio was varied by controlling the degree 
of anneal. One can only conclude that either the scattering due to lattice 
deformation is energy dependent, or some energy dependent coupling 
exists between intrinsic and residual resistivity. 

In essence, the original assumption of independence between these 

components of resistance implies the assumption that the density of 

"impurity"- type scattering centers is so low as to have negligible influ- 
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ence on the overall lattice. Recent investigations may indicate the sev- 
erity of this restriction: Several metals (e. g. , Ag, Au, Cu, Mg) exhibit 
an anomalous resistance minimum, below which resistance increases as 
temperature is decreased (30, 44-46). MacDonald and Pearson (47), 
working with extremely dilute copper alloys, traced this behavior to the 
presence of specific impurities: They reported that as few as 5-10 p. p. m. 
of such elements as Fe, Bi, Pb, Sn, Ga, In, Ge, and Si sufficed to pro- 
duce the phenomenon. In the absence of these elements, no minimum was 
observed. (A later paper (48) describes iron as the primary culprit, 
suggesting that the latter five elements listed serve primarily to replace 
traces of Fe originally tied up in oxides, thus releasing them into solid 
solution. ) 

3. Attempts to "improve the fit" between Bloch's law and experimental 

data such as that proposed by Gruneisen do not alter the asymptotic 

behavior predicted by that law at low temperatures. However, it has been 

5 

found that several metals deviate from the predicted T dependency: for 

2 

example, the resistivity of platinum apparently goes to zero as T (49, 50), 

3 

and that of iron and nickel, as T (51). Perhaps the most startling devia- 
tion is that reported by Mendoza and Thomas (52): Working with copper 
and gold specimens which did not exhibit the minima discussed above 
(which occurs at about 10° K for copper), they found "normal" behavior 
down to about 0. 25° K, below which the resistance of both metals rose 
asymptotically toward the ordinate as temperature was further decreased. 

4. Although Kohler predicted that magnetoresistance should saturate 

at sufficiently high fields, only two metals aluminum and indium (37, 

53) have experimentally displayed such saturation, to date. Further, 

though data from various sources correlate rather well for some metals, ' 
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when plotted on Kohler diagrams, other metals display considerable 
"scatter". 

Lest the trees obscure the forest, it should be noted that the behavior 
of most metals conforms reasonably well to that predicted by the theory, 

» 

especially above the liquid hydrogen range. Even in the "exceptional 
cases", the absolute magnitude of the deviations is small, if only because 
total resistance is small at the temperatures in question. However, from 

a practical viewpoint, a relative error of a factor of two or even, 10% 

might be the difference between success or failure for a particular 

cryogenic magnet. Thus, it is unwise to predict quantitative results, a 
priori, from the theory alone. 

Because of these considerations, the decision was made to conduct 
a careful, detailed survey of the literature regarding the experimental 
resistive behavior of various cryogenic conductors, as a starting point 
for this study. The specific metals investigated were aluminum, copper, 
sodium, and indium: The first three have already been used in the con- 
struction of various cryogenic solenoids, and the last has been proposed 
for use in cryogenic coils, in view of its magnetoresistive behavior. The 
results of this investigation are summarized in the remaining sections of 
this chapter. TO FACILITATE DIRECT COMPARISON BETWEEN CRYO- 
GENIC CONDUCTORS, ALL RESISTIVITY FIGURES IN THIS AND SUB- 
SEQUENT CHAPTERS HAVE BEEN NORMALIZED TO THE RESISTIVITY 
OF "STANDARD" COPPER AT 300° K (WHICH IS TAKEN TO BE 1. 7715 
MICRO-OHM CENTIMETERS), except that data displayed in Figs. 2-5. 

C. "Residual" or "impurity" resistivity : 

As a rough rule of thumb, most metals can be purified sufficiently to 

achieve a residual resistivity approximately four orders of magnitude 
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below their room temperature resistivity. Whether or not such purifi- 
cation is feasible, in a particular case, depends on the cost and complexity 
of the necessary metallurigical processing, as balanced against the con- 
templated application: 

For example, at liquid nitrogen temperatures, residual resistance is 
generally negligible compared to intrinsic resistance, for metals of "com- 
mercial" purity. Thus, the conductors of coils intended for use with this 
coolant need not be of highest purity: Far more importance should be 
assigned to the selection of a conductor with a high Debye temperature, 
to obtain (hopefully) lower intrinsic resistance. (Approximate figures 
for the Debye temperatures of the metals covered by tlis study are: 
indium, 100° K; sodium, 202° K; copper, 333° K; aluminum, 395° K. ) 

By contrast, in the liquid neon range (i. e. , about 30° K) and below, 
residual resistance becomes increasingly important: Ultimately, the 
lowest residual resistance attainable for a given metal generally determines 
the minimum operating power attainable for a coil wound of that metal. 
Obviously, then, the minimum residual resistance attainable, and the ease 
with which it can be attained, become important factors in the evaluation 
of prospective cryogenic conductors. 

In general, metals such as indium and sodium, with low Debye tem- 
peratures, can be readily purified by vacuum distilation: indium, to a 

residual resistance ratio as low as r n /r onn . . 1ri -5 /C/1 v , ,. 

0' 300 = 4. 4 x 10 (54), and sodium, 

to about r o/ r 3oo = 3 x 10 ^ (41, 55, 56). Or, normalizing these figures 

to the resistance standard cited above, the ratio for indium becomes 

i-4 i-4 

r Ql r e = 2 x 10 , and for sodium, r 0 /r g = 8 x 10 . It is worth noting 

that metals with such low Debye temperatures are fully annealed at room • 

temperature: Thus, once residual resistivity has been decreased to a 
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particular level, it will not increase due to "cold work" in the manufacture 
of a coil. 



Aluminum, despite its high Debye temperature, can be readily pur- 
ified by the zone refining process (57-59). This technique has proven a 
relatively economical source of aluminum with residual resistance ratios 
as low as r Q/ r 3Qo = 1* 4 x 10 ( r o' r s = 2. 2 x 10 ). However, both 

Maimoni's data (40) and the Project Sherwood report (20) indicate that the 
resistivity of aluminum is quite sensitive to cold work. If aluminum cryo- 
genic coils are to derive maximum benefit from the use of highly pure 

metal, conductors must be carefully annealed after assembly, if 

possible and structural support provided to ensure that working stress 

in the windings does not exceed the proportional limit. ( It might be noted, 
also, that elimination of impurities tends to reduce the yield strength, 
as well as the resistivity, of most metals. ) 

Copper poses a unique problem: although methods have recently 
been perfected to zone refine this metal (60), its high conductivity com- 
plicates the use of induction heaters in the refining process, thus render- 
ing the technique relatively difficult. To date, most of the ultra-pure 
copper used in recent experiments has been prepared by a complex elect- 
rolytic process developed in 1941 (61). Minimum residual resistance 

-4 -4 

ratios have been in the vicinity of r ()/ r 300 = 5 x 10 , or ^Q/ r s - 4. 8 x 10 

(17, 48, 62). However, though the resistivity of pure copper is apparently 

not so sensitive to cold work as that of aluminum, it is acutely sensitive 

to the presence of Fe in solid solution (as has been discussed above): so 

much so, that actual attainment of the minimum cited is a rare phenomenon. 

(As contrasted with the usual correspondence between decreased impurity 

content and lower residual resistivity, MacDonald and Pearson (47) indicate 

that the resistivity of Fe- contaminated copper can be decreased by reacting 
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the iron contaminant with oxygen, thus raising the total impurity content, 
but converting the iron to less critical iron oxide. ) It would appear advis 
able to avoid even the use of steel dies in drawing pure copper into wire, 
and even then, a copper-wound coil might be a questionable choice for 
low temperature applications. 

The above resistance figures do not, by any means, constitute an 

absolute minimum: Gold, et al. (48) report a specimen of purest natural 

-4 

copper with a resistance ratio of 3. 1 x 10 . Even more significantly. 

Post and Taylor have recently developed a process (63) which, while 
altering neither the purity nor the degree of anneal, lowers the residual 
resistance ratio for sodium from t he level cited above to r()/r^ nn = 1. 5 x 
10 ^ (rQ/r ff - 4. 3 x 10 ^) . It is at least possible that their process may 
be effective for other metals. It might also be noted that the advancing 
state-of-the-art has roughly reduced minimum residual resistance ratios 
by two orders of magnitude since the 1930's. It would seem reasonable to 
assume that at least some further improvement is probable. Again, any 
further reduction in residual resistance will generally permit further 
reduction of the minimum attainable operating cost for cryogenic coils. 

D. "Intrinsic" or "Ideal" resistance: 

Figures 2-5 present a summary of the available experimental data 

on the intrinsic resistance of indium, sodium, copper, and aluminum. 

In processing this data, occasional interpretation has been necessary 

Some of the data were originally presented in the form of curves of "total 1 

resistance, with no specification as to the residual resistance of the 

specimen tested. Such curves have been graphically extrapolated to 

0° K, and the intrinsic resistance computed as discussed earlier. A few 

experimenters have presented curves normalized to their estimate of the 
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Debye temperature resistivity of the specimen tested. In such cases, the 
data have been converted in accord with the Bloch-Gruneisen theory to the 
form presented here: Where the original author specified a value for the 
Debye temperature, that value was used in the conversion. Where no 
value was indicated, the temperatures cited in the preceding section were 
used. 

Since the great bulk of resistivity data for indium, sodium, and copper 
were originally normalized to the 0°C resistivity of the specimen tested, 
the experimental data presented in Figures 2-4 are all normalized to that 
standard. Conversely, published data for aluminum has most frequently 
been normalized to the 300° K resistivity, and is so normalized in Fig. 5. 

For comparison, the Bloch-Gruneisen curve for each metal is plotted 
along with the experimental data. For sodium, copper, and aluminum, 
the coefficients andQ^ used in computing the theoretical curve were 
those determined by Gruneisen (34). For indium, the values were deter- 
mined by the author. 

Not all available data have been included in Figs. 2-5: However, all 
recent works known to the author are included, and the data presented are 
considered fairly representative of the total information available. (In 
fact, it would seem that the trends which appear in the data presented are 
established with sufficient clarity to warrant, perhaps, greater confidence 
in the "overall picture" than in any single work. ) 

As is apparent in Fig. 2, experimental data for indium conform quite 
well to the Bloch-Gruneisen theory down to about 50° K. Below this temp- 
erature, the data lie above the theory: The greatest divergence occurs in 
the range 4° K to 12° K, where some experimental points lie more than 
a factor of 3 above the theoretical prediction. Despite this divergence, the 



experimental curves exhibit approximate parallelism to the theoretical 
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curve: In fact, it would appear though available data are scanty that 

the experimental curves again converge toward the theoretical curve 
below 4° K. Note that, due to its low Debye temperature, the intrinsic 
resistance of indium exceeds the residual resistance figures cited earlier 
down to about 3° -4° K. Considering this, there would be little gained 
by seeking further decrease in residual resistance for this particular 
metal. Indium undergoes a field-free transition to superconductivity at 
3. 41° K. However, since it is a "soft" superconductor, this transition 
would presumably be quenched, in the unlikely event that indium high- 
field coils were designed for operation below that temperature. 

In many respects, the resistive behavior of sodium resembles that 
of indium, as can be seen by comparing Figs. 2 and 3. However, there is 
generally better correspondence between experimental data and the Bloch- 
Gruneisen theory for sodium. Although some divergence is again apparent 
below 50° K, the data of MacDonald and Mendelssohn coincide with the 
theory to 11° K. The greatest divergence is that reported by Woods, at 
about 10° K, and amounts to no more than a factor of 2. Again, the experi- 
mental data appear to reconverge to the Bloch-Gruneisen curve at low 
temperatures. The residual and intrinsic resistivities of pure sodium are 
of roughly equal magnitude around 10°-11° K, in keeping with its higher 
Debye temperature. One would expect that at least some advantage could 
be gained from further reduction of the residual resistivity of sodium, 
since Taylor and Post (26) indicate that "minimum power" sodium coils 
should be operated at about 10° K. 

Fig. 4 indicates remarkably close correlation between experimental 
data and theory for copper, provided that no iron contamination is present: 
in such cases, divergence is generally less than 10%. (The scarcity of data 
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points below 15° K again represents the severity of this restriction!) The 
data of Los and Gerritsen, and of Gerritsen and Linde, represent typical 
behavior in the more normal situation, when iron contamination is present. 
(These curves are labelled "qualitative" behavior, due to the difficulty of 
determining any "constant" value for residual resistance, when total 
resistance rises with decreasing temperature. ) Were it not for the dif- 
ficulty of obtaining iron-free copper, these data would indicate that copper 
might be an ideal choice for low-temperature coils, since its intrinsic 
resistance decreases so rapidly. Actually, the magnetoresistive behavior 
of copper is the least desirable of the four metals considered here, as 
will be shown in a later section of this chapter. These facts, taken together, 
limit the desirability of copper cryogenic coils to the upper end of the 
0°-100° K temperature range under investigation. 

The resistive behavior of aluminum, as depicted in Fig. 5, is unique 

in many respects: First, strong coupling between residual and intrinsic 

resistivity is apparent. The three curves plotted from Maimoni's data 

-4 -4 

represent specimens with r p' r 300 = 476 x 10 , 4. 94 x 10 , and 3. 62 

_ 3 

x 10 , respectively, with direct correspondence between increased 

intrinsic resistivity and increased residual resistivity. This, of course, 
is in direct disagreement with Mathiessen's rule, and is the apparent 
source of the noted strain sensitivity of aluminum conductors. Aside from 
the fact that this behavior results in marked deviation from the Bloch- 
Gruneisen curve (as did Fe- contamination in copper), the nature of that 
deviation is also unique: At the upper end of the temperature span investi- 
gated, aluminum conforms even more closely than did copper to the 
behavior predicted by the theory. There is virtually no deviation between 
experimental data and the theoretical curve. However, at some tem- 
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perature between 65° K and 25° K depending upon the residual resis- 
tance of the specimen the experimental curves shear sharply away 

3 

from the theory, after which they exhibit a T_ dependence, while the 

5 

theoretical curve rapidly approaches the usual T dependence. This 
behavior does not appear to persist for long: Around 10° K, the experi- 
mental curves converge rapidly, and thereafter appear to parallel the 
theoretical curve, though they do not converge to it. (Fenton, et al. , 
reported that their data conformed to the relationship 

P = 5. 68 x 10 -10 + 1. 5 x 10" 13 T 2 * 9 Q-cm 
throughout the range 0 O -30° K. Actually, their data points appear to fall 

increasingly under this line at temperatures below 4° K, and over it at 
temperatures above 25° K. ) Because of this "sidestep", one faces the 
paradoxical situation that predictions from the theory appear to be pre- 
cisely correct, for purest aluminum, down to about 25° K, but more 
than an order of magnitude too low at temperatures much below 25° K. 

It is interesting to speculate that aluminum, above all metals, would 
profit from the further advance in metallurgical technique required to 
further reduce minimum residual resistance levels: Even with the above 
behavior, aluminum has such a high Debye temperature that its intrinsic 
resistance decreases more rapidly with temperature than was the case 
with the other metals covered here. Since the magnetoresistance of 
aluminum apparently saturates at a low value, residual resistance tends 
to dominate the performance of low-temperature aluminum cryogenic 
coils. Further reduction of residual resistance would apparently decrease 
both components of "field-free" resistance, and (if Ar/r Q does indeed sat- 
urate at a constant value) decrease the maximum Ar due to magnetic 
effects, as well. Such a three-fold gain could only result in a direct 
reduction of coil operating cost, again provided that care is taken to ensure 
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that "cold work" does not occur during assembly and operation of the coil. 



One additional conclusion can be drawn from the data plotted in Figs. 
2-5: The Bloch-Gruneisen theory is a useful asymptotic limit for the 
experimentally derived curves. At liquid nitrogen temperatures, data 
points tend to fall on the theoretical curve, or so near it, that the theory 
itself is a perfectly adequate design criterion. But to repeat the obser- 
vation made earlier, the theory should not be used as a design standard 
for cryogenic coils intended for service in the liquid neon range, or 
below. Recourse must be made to experimental data for such applications. 

E. Zero-field resistance standards : 

If the Bloch-Gruneisen theory is not to be used as a design standard 
for cryogenic coils, the time has come to develope a standard. Given 
the above data, this is a relatively simple task, but first, it may be 
profitable to digress to a point of philosophy: 

The primary motivation for use of cryogenic coils for any application 
is the desire to reduce operating costs. Even if a particular cryogenic 
coil is not designed to seek minimum operating cost (for example, geometry 
may be altered from the "minimum cost" optimum to reduce coil size, or 
operating temperature altered, because of preference for a particular 
coolant, etc. ), the primary justification for accepting the greater com- 
plexity and initial cost of cryogenic coils, as compared to "conventional" 
coils, is still the same. 

Bearing this in mind, the "total field-free resistivity" standards 
evolved here will assume that every effort has been made to achieve 
minimum residual resistivity, in each case. As was pointed out above, 
this probably would not be done for liquid-nitrogen- cooled coils: However, 
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intrinsic resistance tends to dominate the performance of such coils to 
such an extent that the error incurred by assuming too low a value for 
residual resistivity is negligible. If, on the other hand, a coil is designed 
for operation at temperatures where residual resistance is a significant 
part of total field-free resistance, the above assumption is strongly justi- 
fied: At such temperatures, total operating cost is dominated by refrig- 
eration costs. For every unit of heat generated in the coil, many units 
of work must be expended in refrigeration. It is highly probable that the 
additional cost and effort required to obtain and maintain minimum residual 
resistance would be repaid many times over the operating life of the coil. 
The importance of this fact cannot be overemphasized! 

For this reason, the residual components of the resistance standards 
presented in Tables 1-4 are the minimum figures previously cited for each 
metal. The intrinsic component corresponds to the Bloch-Gruneisen 
theory, down to the temperature at which the purest specimens tested 
diverge from that theory, and then follows the experimental trend. All 
figures have been normalized to the 1. 7715 micro-ohm centimeter stand- 
ard. The design standards for all four metals are also plotted in Fig. 6, 
for ease of direct comparison. 

The detailed application of these resistance standards to the analysis 
of cryogenic saddle coils is discussed in Chapter VII, and the mean resis- 
tivity design curves resulting from that analysis are presented in Appendix 
B. 
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TABLE 1 



FIELD FREE RESISTIVITY OF INDIUM 



0 


1 


2 


3 


4 


5 


6 


7 


8 


9 


1 . 987 


1 . 987 


2 . 00 


2 . 15 


3 . 52 


6 . 86 


13 . 44 


25 . 2 


44 . 4 


69 . 4 


103 . 3 ‘ 


146 . 2 


197 . 8 


258 


326 


404 


485 


576 


675 


779 


892 


1002 


1129 


1255 


1381 


1507 


1642 


1773 


1913 


2060 


2200 


2350 


2500 


2650 


2790 


2940 


3090 


3240 


3390 


3540 


3690 


3840 


3990 


4140 


4300 


4450 


4600 


4760 


4910 


5060 


5220 


5380 


5540 


5700 


5850 


6010 


6170 


6330 


6490 


6650 


6800 


6960 


7110 


7280 


7440 


7600 


7750 


7810 


8060 


8220 


8380 


8540 


8700 


8860 


9010 


9160 


9320 


9480 


9640 


9800 


9950 


10110 


10260 


10420 


10570 


10730 


10880 


11040 


11190 


11350 


11510 


11660 


11820 


11980 


12140 


12300 


12460 


12620 


12780 


12960 



Values of 10^r/r are tabulated vs. temperature in degrees Kelvin, 
s 

Top temperature scale: units. Left temperature scale: tens. 
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TABLE 2 



FIELD FREE RESISTIVITY OF SODIUM 



0 


1 


2 


3 


4 


5 


6 


7 


8 


9 


4. 33 


4. 33 


4. 33 


4. 33 


4. 35 


4. 40 


4. 51 


4. 73 


5. 11 


5. 75 


6. 74 


8. 15 


10. 30 


13. 23 


17. 10 


22. 3 


29. 1 


37. 2 


47. 9 


59. 8 


75. 9 


94. 8 


115. 1 


140. 0 


165. 7 


195. 9 


230 


268 


308 


352 


399 


450 


504 


559 


619 


680 


743 


810 


880 


953 


1027 


1103 


1182 


1262 


1343 


1424 


1507 


1590 


1673 


1758 


1843 


1930 


2020 


2110 


2200 


2300 


2390 


2490 


2580 


2670 


2760 


2850 


2940 


3040 


3130 


3230 


3330 


3430 


3530 


3630 


3730 


3830 


3930 


4030 


4130 


42 30 


4330 


4430 


4530 


4630 


4730 


4830 


4930 


5030 


5130 


5230 


5330 


5420 


5510 


5610 


5710 


5810 


5910 


6010 


6110 


6210 


6310 


6410 


6510 


6600 



Values of 10^r/r g are tabulated vs. temperature in degrees Kelvin. 
Top temperature scale: units. Left temperature scale: tens. 
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TABLE 4 



FIELD FREE RESISTIVITY OF ALUMINUM 



0 


1 


2 


3 


4 


5 


6 


7 


8 


9 


2. 16 


2. 16 


2. 17 


2. 18 


2. 20 


2. 25 


2. 31 


2. 40 


2. 51 


2. 65 


2. 83 


3. 02 


3. 22 


3. 46 


3. 72 


4. 02 


4. 37 


4. 79 


5. 27 


5. 84 


6. 52 


7. 30 


8. 23 


9. 39 


10. 80 


12. 67 


14. 91 


17. 61 


20. 6 


24. 1 


28. 1 


32. 7 


37. 8 


43. 7 


50. 1 


57. 3 


65. 3 


74. 1 


84. 0 


94. 9 


106. 7 


119. 3 


133. 2 


147. 8 


163. 0 


179. 9 


197. 0 


215 


234 


255 


277 


300 


324 


350 


378 


407 


436 


467 


498 


529 


561 


595 


630 


666 


703 


741 


781 


822 


865 


908 


952 


998 


1044 


1091 


1138 


1185 


1232 


1279 


1326 


1374 


1422 


1470 


1519 


1569 


1620 


1673 


1727 


1781 


1835 


1890 


1946 


2003 


2060 


2117 


2174 


2231 


2288 


2345 


2402 


2459 






4 , 

















Values of lO^r/r are tabulated vs. temperature in degrees Kelvin, 
s 

Top temperature scale: units. Left temperature scale: tens. 
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FIGURE 6 



Field-Free Resistance Standards 




000 li t H 



F. Magnetoresistance: 



Basically, the comments in Section D of this chapter regarding the 
processing of experimental data apply as well to this section. In addition, 
a few of the problems associated with the collection of experimental 
magnetoresistance data bear discussion, prior to presentation of that data. 

Assuming that Kohler's rule is valid, the most convenient manner to 
collect magnetoresistance data is to obtain the purest specimen available, 
immerse it in liquid helium, and subject it to various magnetic fields 
while conducting resistance readings: The purity and low temperature 
will assure a low value for r , and thus a high value of H/r for rela- 
tively low values of H. These data can then be projected to much higher 
field intensities, for coils operating above the temperature of liquid 
helium. 

Unfortunately, aside from the obvious problems of instrumentation 
and measurement, and the difficulties inherent to any work in liquid 
helium, such collection of magnetoresistance data is beset by a number 
of unique problems: 

As zero-field resistivity is quite low, it would be convenient to work 
with a relatively long, thin wire to maintain conveniently large resistance 
values. However, size effects are known to influence magnetoresistance 
readings: The resistance of a fine wire actually decreases in the presence 
of a magnetic field, until the field intensity reaches the point that the 
Larmour radius of a free electron in the field is of the order of the radius 
of the wire. Beyond this point, magnetoresistive behavior is more in 
keeping with the theory. Still, the extent to which size effects may con- 
tinue to influence behavior, even beyond this crossover point, is difficult 
to determine. Unfortunately, the working area available to the experi- 
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menter is obviously limited by the size of both the dewar and the magnet 
providing the test field. Due to these problems, very little data have 
been taken from "large" specimens: that is, specimens approaching the 
size of conductors which might be used in cryogenic coils. In fact, the 
very small size of some specimens tested leads to problems in correlat- 
ing the available data, and questions as to whether such data can be accur- 
ately projected to coil performance. 

Orientation is a second problem: it is known that the magnetoresis- 

tance of a polycrystalline specimen rises more rapidly if the field is 
applied transverse to the specimen, than if the field is applied longitud- 
inally. The data presented here are all transverse-field data, in keeping 
with the anticipated application. Yet this factor also introduces an obvious 
element of uncertainty. 

Finally, the range of the data becomes a problem: Specimens are not 
always so pure and strain-free as might be desired, and attainment of 
sufficiently high values of H/r Q can be a problem. The usual solution is 
to use a pulsed magnetic field, since few magnets exist which are capable 
of producing the necessary field intensities on a steady-state basis. Such 
dynamic readings introduce yet another element of uncertainty: The spe- 
cific heat of most metals is quite low at the temperature in question. 

Pulsed fields induce eddy currents in the specimen; eddy currents result 

in ohmic heating; and the temperature of the specimen and therefore 

its resistivity may rise from this effect. Some experimenters have 

corrected their data for such factors, or have estimated the effect to be 
negligible. Others do not even discuss the problem. One would anticipate 
that such an effect might lead to artificially high values of Ar, particularly 
at high field intensities, but the degree of correction required to compensate 
for this may be a moot point. 

In any event, it is apparent that^the collection of magnetoresistance 
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data is a somewhat more arduous undertaking than the collection of zero- 
field resistance data, so it is perhaps not surprising that the data dis- 
played in Figs. 7-10 exhibit somewhat more scatter than was apparent in 
Figs. 2-5! 

Figures 7 and 10 reveal marked similarity in the magnetoresistive 

behavior of indium and aluminum: For both metals, Ar^/r n apparently 

saturates at H/r = 10^ KG (where r = r/r ), and thus never rises much 
n ns 

above the value of 1. 5. Further, though sources of data for indium are 
relatively limited, it is generally true for both metals that experimental 
data from different sources correlate reasonably well. (Purcell's data, 
as published in NBS Report 8481 (20), are an exception to this: His data, 
on aluminum foils, exhibit significantly higher values of Ar/r Q than have 
been reported elsewhere. Further, the magnetoresistance of Purcell's 
specimens apparently increased with increasing temperature, for a 
given field strength. Due to this latter phenomenon, his data cannot be 
reduced to a single curve by application of Kohler's law, and are not 
reproduced here. The extent to which size effects may be responsible 
for such behavior is not known. ) 

Of the data reproduced in Fig. 10, the only significant deviation 
from the trend cited above is that reported by Cotti, et al. : Their data 
points indicate a possible departure from saturated behavior at high 
field intensities. Other experimenters have generally failed to observe 
such a phenomenon, and Cotti describes any such departure as "spurious". 

II 

In particular, Callaghan whose 3/8 diameter specimen may be the 

largest aluminum specimen for which magnetoresistance figures are 

available tested to virtually the same range, with no apparent deviation 

from saturation. 

Figure 8 presents the experimental data for sodium, and is, in itself, 
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FIGURE 9 



Magnetoresistance of Copper 
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ample evidence of the "bizarre behavior" alluded to earlier. Hopefully, 
the situation may not be quite so chaotic as it would seem, at first. 

White and Woods worked with very fine (less than . 35 mm) rods of 
sodium, and as a result, the magnetoresistance of their specimens was 
found, initially, to decrease with increasing field. The data points 
reproduced here are limited to those which fell above the "crossover 
point" discussed earlier, but again, the extent to which size effects might 
still exert an influence is not known. In any event, the experimental 
curve exhibits a "knee" at about H/ r^ = 300 KG, as though magnetoresis- 
tance had begun to saturate, but then passes through a point of inflection 
and asymptotically approaches a linear dependency at high fields. Mac- 
Donald's data overlap the upper end of this range, and confirm the linear 
dependency at high fields: however, MacDonald's curve is displaced lat- 
erally from that of White and Woods. This may, in part, be attributable 
to the fact that MacDonald indicates his calibration to be "approximate". 

The most recent data reproduced here are from Taylor and Post: 
Their curve results from measurements on a particularly large (sodium 
cast into a 20' length of 0. 5" O. D. x 0. 035" wall stainless tubing) and 
partucularly pure 2 ~~ specimen, wound into a 4-1/2" 

diameter coil. Qualitatively, the behavior of this larger specimen is 
similar to that described above. The curve exhibits a "knee", though a 
much more gradual one, and apparently passes through a point of in- 
flection near the upper end of the range investigated. Whether or not the 
quantitative differences are attributable to size effects is a matter for 
speculation. Regrettably, current data do not extend to sufficiently high 
values of H/r^ to confirm whether or not this last curve would eventually 
approach a linear field dependency. Taylor indicates that additional test- 
ing will extend the data to higher fields. 
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Finally, Fig. 9 summarizes the experimental data for copper. Con- 
trasted with the other metals described here, the magnetoresistance 

2 

curve for copper seems almost simple: It initially exhibits an H depend- 
ence, then falls gradually to a linear dependence. This fact has been con- 
firmed by numerous investigators, and all reports correlate remarkably 
well, except for that of Olsen and Rinderer. Although the magnetoresis- 
tance of copper is relatively low at low fields, it rises so rapidly that 
4 

by H/r^ = 10 KG, copper exhibits greater magnetoresistance than all 
other metals considered. As ultra-low-temperature, ultra-high-purity 
high-field copper coils would exceed this value of H/r n by an appreciable 
margin, one faces the problem discussed above: High-purity copper 
coils are unattractive in the liquid hydrogen range because of magnetoresis- 
tance. Coils of lesser purity are unattractive in that range due to the 
overwhelming probability of iron contamination. Thus, it is advisable 
to restrict the use of copper coils to the liquid nitrogen range and above. 

The design standards derived from Figs. 7-10 are presented in 
Table 5, and plotted in Fig. 11. In order to permit analysis of high- 

field coil performance in the liquid helium range, it has been necessary 

0 

to extrapolate the data to H/r =10 KG. Indium and aluminum have been 

n 

extrapolated at constant Ar; copper has been extrapolated under the 
assumption that Ar varies linearly with field intensity; and sodium has 
been extrapolated under the assumption that magnetoresistance asymptoti- 
cally approaches a linear field dependency. Extrapolated values are indi- 
cated by parentheses in the table, and broken lines in the curves. 
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TABLE 5 



MAGNETORESISTANCE DESIGN STANDARDS 



H/r n (KG) 

1 x 10 1 

2 
3 



x 10 3 



In: Ar/r 


Na: Ar/r 


. 00150 


:: 


. 00270 


-- 


. 00409 





. 00565 




. 00730 




. 00855 


_ _ 


. 00949 




. 01067 




. 01228 


. 00100 


. 0418 


. 00390 


. 0811 


. 00830 


. 1232 


. 01352 


. 1660 


. 01915 


. 212 


. 0270 


. 258 


. 0316 


. 304 


. 0382 


. 349 


. 0447 


. 393 


. 0513 


. 762 


. 1101 


1. 020 


. 1613 


1. 205 


. 206 


1. 340 


. 243 


1. 441 


. 276 


1. 523 


. 306 


1. 589 


. 334 


1. 640 


. 358 



Cu: Ar/r Al: Ar/r 



— 


. 00110 





. 00197 




. 00309 




. 00448 


. 00140 


. 00606 


. 00182 


. 00794 


. 00229 


. 01008 


. 00284 


. 01234 


. 01093 


. 0448 


. 0230 


. 0840 


. 0378 


. 1268 


. 0545 


. 1701 


. 0716 


. 214 


. 0894 


. 256 


. 1080 


. 298 


. 1273 


. 339 


. 1470 


. 563 


. 364 


. 676 


. 599 


. 872 


. 838 


1. 008 


1. Ill 


1. 112 


1. 342 


1. 187 


1. 591 


1. 249 


1. 844 


1. 295 


2. 11 


1. 338 



Continued on next page 
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TABLE 5 



(CONTINUED) 



H/r (KG) 
n 

4 

1 x 10* 

2 
3 



x 10 5 



In: Ar/ r 


Na: Ar/r 


Cu: Ar/r 


Al: Ar/r 


1. 668 


. 381 


2. 36 


1. 369 


1. 870 


. 552 


4. 94 


1. 489 


1. 906 


. 704 


(7. 41) 


1. 508 


(1. 906) 


. 842 


(9. 86) 


1. 513 


(1. 906) 


. 969 


(12. 30) 


(1. 513) 


(1. 906) 


(1. 097) 


(14. 77) 


(1. 513) 


(1. 906) 


(1. 217) 


(17. 30) 


(1. 513) 


(1. 906) 


(1. 330) 


(19. 78) 


(1. 513) 


(1. 906) 


(1. 440) 


(22. 2) 


(1. 513) 


(1. 906) 


(1. 552) 


(24. 7) 


(1. 513) 


(1. 906) 


(2. 54) 


(49. 5) 


(1. 513) 


(1. 906) 


(3. 42) 


(74. 4) 


(1. 513) 


(1. 906) 


(4. 26) 


(99. 6) 


(1. 513) 


(1. 906) 


(5. 06) 


(124. 9) 


(1. 513) 


(1. 906) 


(5. 86) 


(150. 0) 


(1. 513) 


(1. 906) 


(6. 67) 


(174. 7) 


(1. 513) 


(1. 906) 


(7. 47) 


(199. 6) 


(1. 513) 


(1. 906) 


(8. 32) 


(225) 


(1. 513) 


r = 


r/r 






n 


s 






r = 


1. 7715 x 10" 6 


ohm -centimeters 





53 



AR/R 



FIGURE 1 1 



Magnetoresistance Standards 




H(R/R) (KG) 
s 

54 



AR/R 



III. STRUCTURAL DESIGN CONSIDERATIONS 



A. Basic concepts : 

The non-conducting structural components of any air- core electro- 
magnet must be regarded as a necessary evil: While necessary to support 
the coil against the various stresses developed, such components con- 
tribute nothing to the desired objective of generating a specified field 
throughout a specified region. Suffice it to say that if coils could be made 
to operate without such elements, they would be. Since this is generally 
impossible, this chapter will be devoted to the development of a model 
with which the minimum structural requirements of a cryogenic saddle 
coil may be estimated, as a function of the basic variables defining the 
coil. 

Note that the above does not imply that the model used provides 
answers which can be directly applied to "real" coils. For one thing, 
the minimum structural fraction achievable in real coils is a function of 
such things as the ingenuity of the designer, the need to provide for ease 
of construction, etc. For another, a model which determines the min- 
imum structure "needed" for support may quite reasonably predict that 
no structure is necessary, in a particular instance. This, for example, 
would be true for small, "jelly-roll" cryogenic solenoids, in which the 
configuration permits the conductors to be essentially self-supporting. 
However, even in small "jelly-roll" coils, it is convenient to use non- 
conducting elements to ensure centering, provide for coolant flow, etc. 
Thus, the structural volume fraction of a "real" coil may not be large, 
but will not be zero. 

Nevertheless, the model proposed here provides answers which are 

useful in the preliminary design phase, which is consistent with the spirit 
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of this work. No attempt will be made to evaluate those variables which 
are primarily a function of the discretion of the designer. Corrections 
for such items can best be injected during the detailed design phase, after 
preliminary analysis has broadly defined the form of the desired coil. 

In proceeding with this analysis, it will be convenient to discriminate 
between two classes of structural elements: 

Internal structure is here defined as those structural components 
contained entirely within the coil dewar (and therefore, within the coil 
windings). The primary influence of such structure on the design is to 
reduce the packing factor attainable, and thus increase the conductor 
current required to obtain a given effective current density. (Note that 
equations (1), (3), and (5) essentially define the form of a coil in terms 
of effective current density. ) 

External structure is defined as those structural components which 
are external to the coil windings, including the dewar and those structural 
elements which penetrate the dewar. External structure has no influence 
on the effective current density, but does decrease the available working 
area, add to coil weight, and contribute to heat leak. 

B. I nternal structure: 

There are two significant differences between the internal structure 
of a saddle coil and that of a solenoid, such as those described earlier: 

First, the conductors of a tape-wound or "jelly roll" solenoid have 
structural continuity in the plane of main stress, while those of a saddle 
coil do not. The structural strength of the conductors of a tape-wound 
coil can directly support part or all of the stress load. By contrast, the 
conductors of a saddle coil run axially, while the major stresses occur 
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in the transverse plane. The latter situation is clearly less favorable. 

Second, the magnetic stresses developed in a high field solenoid tend 
to pull the coil apart. Thus, it is relatively easy to design internal 
structural elements so that all structural components are in tension. The 
normalized magnetic stress distribution for various saddle coils is shown 
in Figs. 12-16. (These figures also indicate coil profile, as a function of 
the form factor defined in eqn. (6). ) As can readily be seen, it is impos- 
sible to support such a stress distribution by components entirely in ten- 
sion. The point, here, is that structural components bearing compressive 

loads may be subject to elastic instabilities i. e. , buckling at stress 

levels well under the yield stress of the material used. Allowance must 
be made for this in coil design. 

Two other considerations apply to the internal structure of cryogenic 
coils, in general. Mention has been made of the necessity to minimize 
residual resistance for ultra-low -temperature coils. This might neces- 
sitate designing to allow for a final anneal after the coil has been assembled. 

(Or, if the coil is to be built up from modular units that is, pre-formed 

"bundles" of conductors designing for anneal of such units, after shap- 

ing. ) Such goals probably preclude the use of any structural materials 
which depend on elaborate heat treatment to develope structural strength. 
However, the designer has one consolation if "as annealed" metals must 
be used for internal structure. Engineering metals which do not become 
brittle at cryogenic temperatures exhibit an appreciable increase in 
strength, over room-temperature characteristics. (See McClintock and 
Hauser (89) for a discussion of factors influencing mechanical properties. ) 
Since internal structure will be maintained entirely at cryogenic temper- 
atures while the coil is operating and magnetic stresses do not exist 

while the coil is not operating advantage can be taken of the increased 
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low-temperature strength of structural metals. (For example, Swenson 
(90) found that the compressive yield strength of Type 304 stainless steel 
increased from about 75, 000 p. s. i. at room temperature to over 110, 000 
p. s. i. at 20° K. For additional information, see also (91). ) 
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C. Superposition model for internal structure: 



As is apparent in Figs. 12-16, the stress distribution in a high-field 
saddle coil is rather complex. It is probable that detailed structural 
design will be a difficult, iterative analysis, with the final outcome 
strongly dependent upon the arrangement scheme initially assumed. In 
order to reduce this situation to one which is sufficiently tractable for 
preliminary analysis, a superposition model will be used. Basically, the 
illustrated stress pattern results from two separate force fields: that due 
to the interaction between the current and magnetic field of a single coil 
lobe, and that due to the interaction of the current of one lobe with the 
field of the other. Each of these force fields will be treated separately, 
and the optimum structural arrangement to support each field will be 
determined. The structural volume fraction required for the combined 
fields will be estimated as the sum of the structural volume fractions so 
determined. 

As will be seen, this is equivalent to assuming that the support 
structure of the coil is fabricated from a series of interlocking segments 
of circular shells. This probably would not be feasible, for a "real" coil, 
due to the obvious problems of assembly. However, the circular shell 
model does indicate those regions of the coil which most need structural 
support, and provides a tentative arrangement which could be successively 
modified in an orderly fashion. 

Specifically, application of the superposition technique to this problem 
involves the following approximations: 

1. Interactions between the two force fields are ignored. (If one were 

to draw an "x" through the center of either cylinder, the two force fields 

would be found, generally, to oppose each other on that side of the "X" away 

from the channel. Thus, less structure would actually be required in that 
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region than the model would predict. ) 

2. Interactions between support shells are ignored. (The shells 
required to support "self field" are subjected to compressive loading, 
and must be designed in accord with the appropriate elastic stability 
criteria. Actually, the interlocking shells resulting from the "adjacent 
field" analysis would provide some additional stiffening, not considered 
in the "self field" analysis. On the other hand, if either set of shells 
were to shift position slightly under load, part of that load would be trans- 
ferred to the other set. Such load transferral is also neglected. ) 

Basically, predictions from such a model should be reasonable, but 
conservative. At least, the model provides a systematic means to esti- 
mate the additional structure required due to the inability to avoid com- 
pressive loading. Further, it is well to be conservative at this point. 

Deviations from the idealized circular shell model in the interest of 

ease of assembly will probably necessitate a somewhat higher struc- 

tural volume fraction than would result from an "optimum support" arrange- 
ment. 

D. Support of "self field ": 

The force distribution in a circular cylinder, due to interaction of 
uniform, axial current flow with the magnetic field generated by that 
current flow, is directed radially inward and has a magnitude given by: 

f m = magnetic force density = ^ (10) 

To return to our original hypothesis, where the use of structure can 
be avoided, it will be. If the cylinder is comprised of many individual 
conductors, those near the center of the cylinder are in such a region. 

Force densities are low, and directed so as to hold the coil together. It 
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is reasonable to assume that the wires in any circular shell would tend to 
"jam together", as do the keyed stones of a circular arch. Certainly, the 
conductors could be shaped to accentuate this tendency. Stability would 
not be a problem, since this region is essentially solidly filled with con- 
ductors. To what radius does this region extend? 

Consider a thin shell of conductors at radius r. If we assume this 
shell to be "keyed together", it is essentially self-supporting. The mag- 
netic force distributed within this shell is equivalent to an external pres- 
sure, directed inward, of magnitude 



P o j2 r dr 



Neglecting stability, this is simply the classic "hoop stress" problem. 

If we desire to limit the internal stress in the conductors to a maximum 

r P, 



level , we have: 






^J 2 r 2 dr 



t 



2 dr 



and therefore 

r 

max 

That is, the individual conductors must not be "keyed together" at 

r>r , or the resultant hoop stress will exceed the specified stress 
max 

limit of the conductor. If the coil must be designed to prevent self support 
in shells at greater radii, the conductors in such shells will be forced in- 
ward by the magnetic pressure, due to lack of structural continuity. The 

resultant stress must not be transferred to the shell at r = r , since 

max 

the conductors of that shell have already been "stressed to the limit". 
Therefore, the first structural support shell must be located at r = r 
Assuming the conductors located outside the first support shell all 
transfer their load inward, those located immediately adjacent to the sup- 
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port shell will be subjected to the greatest load, in the form of simple 
radial compression. Again, this must not exceed the specified limit, 

. The position of the second support shell can then be determined by 
integrating eqn. (10) radially, subject to this condition. 



. AF 



In a slab of width AO, 

<A rj2 > (r dr de > /V 2 , s 3> , „ 

= (r 2 -rj ) A© 



AF = 




2 6 
The "equivalent pressure" on the structural shell 



is thus: 



l] j2 ii j2 

P e = AF/ AA =— — (r 2 3 -r 1 3 )A©/r 1 A© = 



■ 3 3 

r 2 " r l 



( 12 ) 



And the radius of the second structural shell must be: 

6 O' ^ 1/3 



r 2 r l 



1 + 



1 1 t 2 2 
r 



i j 



(13) 



Obviously, this can be repeated until the outer radius of the cylinder 
is reached. Eqn. (13) is used to determine the radius of each successive 
shell, since it is known that p g = 0^ in eqn. (12). For the last shell, eqn. 
(12) is solved for the effective pressure on the shell, by setting r 9 = r , 
the outer radius of the cylinder. 

The above analysis of "effective pressure" has ignored the thickness 
of the support shells. The reason is this: Consider a slab of height h, 
1 - mounted on a zero-thickness shell at y = 0, and carrying 



h a uniform current in the -z direction. In the presence 

of a uniform, x-directed magnetic field, B x , the slab 

will exert a pressure p g = B x< J e h on the support. Now let the support 
have non-zero thickness, so that it extends from y = 0 to y = t. If we 
require the same total current in the slab, the current density must be- 
come J = J e /A > where A = (h-t)h. The pressure exerted on the support 
is now given by p p = B x (J e /A)(hA) = ®x^e^' aS ^ e ^ ore ‘ An analogous 
situation exists in the "self field" problem above, since the magnetic field 
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at a given radius in such a composite cylindrical conductor is primarily 
determined by the "effective" current flow within that radius. Thus, to a 
first approximation, the location of the "next" shell is not influenced by 
the non-zero thickness of the "present" shell, and the appropriate current 
density to use in eqns. 10-13 is J e , the effective current density of the 
coil. Obviously, field perturbations due to the presence of non-conducting 
structure become more severe, and therefore the accuracy of this model 
decreases, asA-*-0. (In fact, when the stability criteria are applied to 
determine shell thickness, it is quite possible to compute a "structural 
volume fraction" greater than unity. ) However, if the analysis indicates 
that support structure alone would result in a packing factor less than, 
say, 0. 25, the coil is obviously not of interest. 

Once the position of each shell has been determined, its required 

thickness must be found. Windenburg (92) summarized, in 1947, the 

equations required for this process. They are: 

/ 2 t 



L<L 



crit \ 



L>L 



P = — 0 - 
1.05 y D 

2. 42E 

p = d -(f ) 314 

2E 



(t/D) 



5/2 



P = 



crit < 



L/D-. 045(t/D 
(t/D) 3 



1/2 



1. 27 \1 -fJC) 

p = 2. 34 O — -1. 06( 0 /E) 1/2 

y D y 



(14) 

(15) 

(16) 
(17) 



Equations (14) and (15) apply to "short" cylinders, while (16) and (17) 
apply to "long" cylinders. Essentially, if a cylinder is very short, elastic 

instability is no problem, and (14) which is a slightly more conservative 

version of the hoop stress equation applies. As the cylinder is length- 

ened, a simple, two-lobe instability collapse is possible at stress levels 
below yield. The "critical pressure" in this region is a function of the 
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elastic modulus, t/D, and L/D, as shown in (15). At still greater lengths, 
the ends are so remote that further removal has no significant influence 
on stability. At this point, (16) applies. However, as (16) was derived 
from elastic stability considerations, it does not recognize that a pro- 
portional limit exists. Thus, collapse can occur from yielding, vice 
instability, at pressures lower than predicted from (16). The appropriate 
equation to predict this mode of failure is (17). 

In using these equations to determine shell thickness, one must 
obviously know p, t/D, L/D, the design stress limit specified for the 
structural material (which would presumably be an order of magnitude 
higher than that established for the conductor), and the elastic modulus 
and Poisson's ratio for the structural material. For "short" cylinders, 
one would solve both (14) and (15). The larger (t/D) so obtained would be 
correct. For "long" cylinders, (16) and (17) would be applied in similar 
fashion. To aid in distinguishing "short" from "long", Windenburg defines 



If X= 5, (14) and (15) apply. If X = 20, (16) and (17) apply. If 5 < 

X <20, try both (14) and (15), and (16) and (17), as above, then keep the 
smaller of the two resultant values for (t/D). 

One problem arises, regarding the crossover from (16) to (17) for 



the parameter: 





P 



e 



(A) 



high-yield- strength metals. For a given 
metal (and a particular L and D), the simal- 



taneous solution of these two equations in- 
volves a cubic in (t/D). Such a solution can 



assume one of three forms, illustrated in 



Fig. 17. In case (A), there are three real, 
unequal roots. In (B), there are still three 
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Fig. 17 








i 






real roots, but two are identical. In (C), there is only one real root, and 
two imaginary roots (i. e. , no physical solution). Equations (16) and (17) 
were originally derived for mild steel. Equation (17) is an attempt to 
express case (B) in terms of the physical parameters of any metal, and 
is successful for mild steels. The physical solution in this case follows 

(16) to the point of tangency, then (17) at higher pressures. Unfortunately, 
for yield strengths in excess of about 50, 000 p. s. i. , substitution of the 
appropriate parameters into (16) and (17) results in case (C): That is, 

(17) no longer provides the equation of a straight line, tangent to (16) at 
the desired crossover point. 

To circumvent this situation, it is necessary to re-examine the phy- 
sical reasons underlying the necessity for "crossover". Eqn. (16) becomes 
invalid at high pressures, for a particular goemetry, because "real" 
metals reach their proportional limit and thereafter deviate from Hooke's 
Law. Thus, (16) predicts too low a value for (t/D), since it assumes a 
constant elastic modulus. Eqn. (17) attempts to correct for this by 
assuming that shell thickness must be increased in direct proportion to 
pressure, beyond a specified point. This specified crossover point would 
presumably lie somewhere between the proportional limit and the yield 
point (2% offset), depending on the characteristics of the metal used. 

A conservative approach will be used here. Following the solution 
of (16), the value of (t/D) so obtained, along with the known value of P , 
will be inserted in the hoop stress equation to determine whether or not 
the resultant stress in the shell exceeds the proportional limit of the 
metal in question. If so, the effective pressure corresponding to the 
proportional limit for the given shell will be determined, and the slope 
of eqn. (16) will then be found at that pressure. Assuming this to be the 
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desired crossover point, these data will be used to generate the equation 
of a straight line, tangent to (16) at the pressure corresponding to the 
proportional limit. The projection of this straight line will then be used 
to determine a corrected value for (t/D) at the effective pressure. 

Finally, after the radius and thickness of each shell have been deter- 
mined, the total volume of these "compressive" shells can be determined, 
and divided by the volume of the composite cylindrical conductor. It will 
be assumed here that the compressive volume fraction of the "crescent 
moon" of a saddle coil lobe is approximately equal to the compressive 
volume fraction of a cylinder of the same radius, carrying the same 
effective current. 



E. Support of "adjacent field": 

The force distribution in either coil lobe, due to interaction of uni- 
form, axial current flow with the magnetic field generated by a similar, 
but oppositely directed, current flow in the adjacent lobe, is directed 
radially outward from the center of the adjacent lobe and has a magnitude 
given by: 



m 



., 2 t 2 
/Vo J 

2r 



where r Q is the outer radius of the adjacent lobe. A natural structural 
arrangement to support such a force distribution consists of a series of 
circular arcs, concentric with the adjacent lobe. The problem here is 
similar to that preceding. The position, thickness, and in this case, 
length of each such arc must be determined. Since this set of structural 
shells is in tension, no stability problems exist. 

To determine the location of the first tensile shell, let us refer to a 
coordinate system centered on the "adjacent" lobe, henceforth called the 
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left lobe. Again considering the lack of 
structural continuity in the transverse plane, 
since the magnetic force considered here 
tends to "pull apart" the individual conduct- 
ors, it will be assumed that the load will be 
transferred radially outward from the center 
of the coordinate system. In the most severely stressed conductors, 
immediately adjacent to the structural shell, we have: 

r - 2 j 2 f( 

o e J 




P = dF/dA = / 2r©f dr/2r 1 © = 
e J ml 

r 



/v 



© dr 



(18) 



2rj© 



/ ,^22 2 
(x - 2a) + y = r Q , 



where 2© is the angle subtended by the circle bounding the right lobe at 
a given radius. Referred to the indicated coordinate system, this circle 
is: 

where 2a is the distance between lobe centers. Further, the equation 

of a circular arc, centered on the origin, of radius r, is: 

2 . 2 2 
x + y = r 

Thus, the x-coordinate of the intercept of the circle bounding the 
right lobe, and above circular arc, is: 



x = 



2 , . 2 2 
r + 4a - r 
o 

4a 



and therefore 




(19) 



( 20 ) 

The above is a bit messy for explicit solution, but can be programmed 
readily for numerical solution. Conceptually, the procedure is the same 
as before. Shells are positioned in accord with the requirement that p g = 
(3 at each shell. Shell thickness is determined from the requirement 



that : 
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where ^ is the design stress limit for the structural material. The 
"tensile" volume fraction then becomes: 



T. V. F. 



E2( C r/o d )r i 2 e i 

' » 



( 22 ) 



where A is the cross sectional area of either lobe, 
c 

Essentially, the above shells transfer the magnetic stresses dis- 
tributed within the windings to the outer edge of the coil: i. e. , the ends 
of such shells must be "anchored" by structure bounding the outer sur- 
face of the lobe. Further, as there is a net repulsive force tending to 
separate the two lobes, provision must be made to tie the bounding 
structure of the two lobes together. However, such structure does not 
penetrate the windings, and therefore has no effect on the packing factor 
of the coil. It is therefore treated as "external structure", as defined 
above. 

It is interesting to note the role of (J and 0^, in determining the 
structural volume fraction predicted by this model. Basically, the num- 
ber of structural shells both tensile and compressive is determined 

by (S Q > and the thickness of each shell by C^. 

The details of the actual computational procedure followed in applying 
this model to the coils investigated are provided in Chapter VII, and the 
results of those calculations, in Appendix C. 



F. External structure : 

The external structure associated with a "real" MHD coil is a suf- 
ficiently strong function of the intended application as to preclude quanti- 
tative discussion in a study of this nature. Nevertheless, a few points 
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bear mentioning, qualitatively: 

The phrase "external structure" has been defined to include the 
structural elements bounding each coil lobe, as well as those members 
connecting the two lobes, the dewar required to maintain the cryogenic 
environment, and the structural members which support the assembled 
coil, both against gravity, and against any drag or thrust generated by 
the MHD channel. 

Most probably, the structure bounding and connecting the coil lobes 
would form the inner vacuum wall of the dewar. Again, advantage could 
be taken of the increased low-temperature strength of materials. 

For coils designed for operation below liquid nitrogen temperatures, 
it is likely that the dewar itself would be of "two-stage" construction. 

That is, everything within the "cold" side of the dewar would be main- 
tained at the coil operating temperature, but a "middle layer" of the de- 
war would interpose a "sandwich" of liquid nitrogen between the cold side 
and ambient temperature. This would reduce heat leak by conduction 
and radiation through the dewar, but far more importantly, allow greatly 
reduced heat leak by conduction along structural elements and electrical 
leads piercing the dewar, into the "cold box". Most of the heat leak 
along such elements could be removed, with reduced refrigeration cost, 
from the liquid nitrogen layer. 

Most probably, the vacuum layers of the dewar would contain either 
a multi-layer "superinsulation" (93-97), or one of the somewhat less 
effective but more readily applied opacified powders (96-98), to minimize 
heat leak. Such construction would ensure that heat leak through the dewar 
would be negligible compared to other losses. (A reasonable estimate, 
given such construction, would be . 05 watts per square meter. ) 
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' From a strictly structural viewpoint, arrangements would presumably 
be made to ensure that heat leak along electrical leads was also negligible. 
McFee (99) indicates a "passive" method of minimizing such heat leak by 
control of lead geometry. "Active" methods light involve using boil-off 
vapor from the cold box to refrigerate the leads. However, for steady- 
state, closed-loop operation, it would be necessary to provide for recon- 
densation of such vapor, and its return to the cold box. 

Probably, the most efficient geometry for support structure which 
pierces the dewar would be the so-called "bicycle-spoke" arrangement, 
in which the coil is supported by a network of long, thin spokes. Careful 
design of such supports will ensure maximum path length from ambient 
to cryogenic temperatures, and therefore minimum heat leak, along 
each support member. Obviously, this is not the most convenient arrange- 
ment from the viewpoint of ease of construction. In small cryogenic 
2 

coils, I R losses dominate so strongly that one might prefer to utilize 
fewer, and more easily installed supports, even at the cost of added heat 
leak. In larger coils, structural leak may be more important. Stekly 
and Woodson (100) provide a means by which the magnitude of such heat 
leak may be readily estimated: Let F be the total force to be supported 
in a given direction (weight, thrust, etc. ). Then, 

A = F/ (f 
s s 

where 0” is the design stress in the support material (note that since 

supports penetrate the dewar, room temperature figures must be used), 

and A is the required effective cross section of support structure. Since 

the rate of heat leak along such structure is given by: 

A 

q = k — - AT, 
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where k is the mean thermal conductivity (see (101) or (102) for typical 



values), L the mean length, and AT the temperature drop across the 
6 

structural elements, we have: 



A first estimate of coil weight is simply F « m c gV, where m £ is the 
mass density of the conductor, and V the total coil volume. Similarly, 
if the structure supporting the MHD channel must penetrate the dewar, 
the associated heat leak can be estimated by setting F =« (p^ - p Q )A ch 
For a matched load, Pg en ~ Fv/2, where v is the mean flow velocity in 
the channel, and an equivalent estimate is F* 2P /v. In any event. 



such relations can be used to estimate the importance of support heat 
leak, and therefore, the importance of optimizing support geometry. 

Since the MHD channel operates at or considerably above room temp- 
erature, the "working area" of a saddle coil is not fully available for 
the intended application. Some of it must be filled by the type of dewar 
and support structure described above. Again, "how much" is a strong 
function of the particular application intended, and is therefore not 
included in any of the curves generated in this study. To compensate for 
this factor, one can enter the coil geometry curves of Appendix A with 
a somewhat greater A ^ than is actually required for the anticipated appli- 
cation. 




(23) 



gen 



76 



IV. REFRIGERATION 



A. Efficiency: 

As was mentioned earlier, the cost of refrigeration is the dominant 
operating cost of ultra-low -temperature cryogenic coils. Such being the 
case, the time has come to discuss refrigerator efficiencies, before pro- 
ceeding further. 

Brief reference to the Carnot cycle will depict the magnitude of the 

problem. In order to maintain a liquid nitrogen environment that is, 

to maintain a refrigerated area at 77. 4° K, given an ambient temperature 

of 300° K an ideal, Carnot refrigerator would require 2. 9 watts of 

input power for every watt of heat removed from the refrigerated area. 

To maintain a liquid helium environment (i. e. , 4. 2° K) under the same 
conditions would require 70. 4 watts of input power for every watt of 
refrigeration provided. As is hardly surprising, considering the working 
environment, "real" cryogenic refrigerators require appreciably more 
power than these "ideal" figures. 

Actually, the spread between "real" and "ideal" performance has 
been reduced significantly, within the last decade. Scott (62) gives a 
description of the various components used in cryogenic refrigerators, 
as well as a discussion of many of the factors limiting refrigerator per- 
formance. Suffice it to say that though the design of a refrigerator may 
be simple in theory, the design of a cryogenic refrigerator is extremely 
complex, in practice. Thus, improved performance has generally been 

attributable to improved technique i. e. , "practical engineering" 

rather than any "breakthrough" in basic concepts. Given the law of 
supply and demand, the effort expended in this field has increased in pro- 
portion to the growth of industrial application of cryogenic technology. 
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Hopefully, further improvement in the performance of cryogenic refrig- 
erators should be expected, under these circumstances. 

Fortunately, it is unnecessary to become intimately familiar with 
the design of cryogenic refrigerators in order to appreciate the two factors 
which most directly influence their relative efficiency (that is, their per- 
formance as compared to the Carnot cycle efficiency). 

First, as with many other mechanical devices, small cryogenic 
refrigerators ('- 1 watt) tend to be considerably less efficient than com- 
parable large refrigerators (over 100 watts). Of course, the degree of 
improvement with increased size eventually "tapers off". Bernert and 
Stekly (103) present curves of relative efficiency vs. size which indicate 
that relative efficiency "saturates" at a power input of approximately 1000 
KW, more or less independently of "cold side" temperature. The signifi- 
ance of this point is that the refrigeration requirements of the coils being 
analyzed here are sufficiently great that all require "large" refrigerators, 
by such a standard. 

Second, relative efficiency is a function of refrigeration temperature. 
The engineering problems which must be solved to maintain continuous 
operation at 4. 2° K are obviously more stringent than those confronted in 
designing for operation at 77. 4° K. Thus, not surprisingly, the relative 
efficiencies which have been attained in helium refrigerators are lower 
than those which have been attained in nitrogen refrigerators. To be 
specific. Fig. 18 has been prepared as a working estimate of the input 
power required, as a function of temperature, for "large", closed loop, 
continuous-operation cryogenic refrigerators. The "actual" curve is 
based on an estimate of the present state-of-the-art, with no allowance 
for anticipated future improvements. The Carnot cycle theoretical limit, 
as well as a few specific data points, are plotted for comparison. 
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FIGURE 18 
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B. Other factors: 



The steady-state refrigeration requirements of the saddle coils 
treated here generally exceed the refrigeration requirements of known 
cryogenic applications. Again, given the law of supply and demand, refrig- 
erators of sufficient capacity to provide for these needs do not fall into 
the category of "on-the- shelf hardware". While it is felt that construction 
of such refrigeration systems is within the capability of existing tech- 
nology, manufacturers are understandably discrete, regarding the speci- 
fications for systems which exist only "on paper". Further, extrapolation 
of size, weight, and cost figures from the region of current practice 
would appear to involve a greater degree of uncertainty than was the case 
with the efficiency figures quoted above. Regrettably, little concrete 
information regarding such factors is available in the open literature. 

In any event, Bernert and Stekly (103) indicate that total refrigeration 
system mass might be estimated from the relationship: 



n 94 

M =15. 7W t u ‘^ 
s 1 

(where is system mass in Kg and W T refrigerator input power, in KW), 
for systems with an input power requirement in excess of 10 KW. In a 
separate article (106), Bernert indicates that capital costs for 4. 2° K 
systems may be estimated from: 



33W r 0,52 = C ^ 20W 0,66 

K S K 

(where) C g is system cost in thousands of dollars, and Wpj is refrigeration 
capacity at 4. 2° K, in watts), for systems with a cooling capacity in excess 
of about 50 watts. Presumably, capital cost for a given cooling capacity 
should decrease as a function of temperature, probably even more rapidly 
than does system mass, considering the reduced complexity of the design 
problem. Certainly, one could scale the above cost estimate in proportion 



80 



to total system mass, to obtain a crude "first estimate" of the capital 
cost of systems operating above 4. 2° K. However, for specific information 
regarding this factor,, as well as the other factors mentioned in this sec- 
tion, one would be advised to turn to the various manufacturers' repre- 
sentatives. 
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V. HEAT TRANSFER 



A. Basic concepts: 

For various reasons, it is convenient to divide the treatment of heat 
removal into two separate areas: Heat transfer, discussed in this chapter, 
deals with the flux of energy from the windings, into the coolant, across 
a fluid-solid interface. (Here, two factors are of particular importance 
to the analysis of cryogenic coils: the amount of heat transfer surface 
required, and the relationship between rate of heat generation, and temp- 
erature rise between bulk coolant, and conductor. ) Coolant flow , discussed 
in the next chapter, presupposes a given heat flux into the fluid stream, 
and explores the influence of this flux on the bulk flow properties of the 
coolant. (Items of interest will be the bulk temperature profile of the 
coolant, the relationship between inlet and outlet pressures, and packing 
factor, etc. ) 

It is also convenient to categorize the overall heat removal process 
as "single-phase", or "two-phase". The former term implies that the 
coolant is a liquid or a vapor, and that energy is removed by sensibly 
heating that liquid or vapor. The latter term implies that the coolant is 
a saturated liquid, and that energy is removed in the form of latent heat, 
by converting a portion of that liquid to vapor. 

The analysis of single-phase heat removal is an inherently simpler 

task than that of two-phase heat removal. For purposes of preliminary 

analysis, one can characterize a single-phase coolant in terms of "mean" 

properties: Particularly, by a mean heat transfer coefficient and a mean 

specific heat. Since energy is removed in the form of sensible heat, the 

latter, coupled with the bulk flow rate, determines the mean temperature 
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rise in the bulk fluid. For the present application, there is an obyious 
interest in keeping this rise as small as practical. Suffice it to say, that 
coupled with the generally low visosities of cryogenic coolants, this leads 
to sufficiently high bulk flow rates as to ensure fully developed turbulent 
flow of the coolant. One can thus determine the mean heat transfer coef- 
ficient experimentally, or utilize an appropriate predictive correlation, 
since such correlations are reasonably accurate under these conditions. 
(Stevenson, for example, based his design of a gaseous-helium- cooled 
solenoid (107) on the McAdams correlation: 

hD/k = 0. 023P °* 4 R °* 8 , (24) 

re’ 

where h is the heat transfer coefficent; D, coolant tube diameter; k, 
thermal conductivity of the coolant; P^, Prandtl number of the flow, and 
R , Reynolds number of the flow. See Rohsenow and Choi (108) for a 
discussion of this and other such correlations. ) Given the mean heat 
transfer coefficient, the information required at this point can readily 
be determined. While it is apparent that considerable search may be 
necessary to determine the "optimum" operating pressure, etc. , the 
procedure is conceptually straightforward. 

By contrast, the analysis of two-phase heat removal is not straight- 
forward: in general, predictive correlations are not sufficiently reliable 
as to permit accurate computation of the heat transfer coefficient. In 

fact, it would seem prudent if not mandatory to construct a precise 

model of any proposed heat transfer scheme, and obtain experimental 
confirmation of any theoretical predictions, before entering the final 
design stage. This would be particularly advisable in the design of large 
cryogenic saddle coils, as the size and complexity of the over-all system 



would render any significant error in the anticipated heat transfer 
coefficient a most costly mistake. 

The above notwithstanding, most of the cryogenic coils which have 
been built, to date, have utilized some form of two-phase cooling. For 
this reason, the rest of this chapter, and most of the next, will be devoted 
to a general discussion of that problem. 

B. Pool boiling charac teristics : 

In a physical sense, the simplest form of two-phase heat transfer is 
that commonly called "pool boiling", in which a heated specimen is immer- 
sed in a large pool of liquid, and cooled by natural- convection boiling, 
with geometries sufficiently "open" to ensure negligible interference with 
convective flow. (These conditions can be approximated with small, 
"jelly-roll" solenoids, which are often cooled by simple immersion. ) 

The solid line in Fig. 19 is a qualitative representation of the sort 
of heat transfer curve which one might typically observe under such 
circumstances. To briefly discuss the regions of this curve: 

Region I represents non-boiling convective cooling, and is character- 
ized by low heat flux, and low heat transfer coefficient. (Note, however, 

that the "AT" referred to in Fig. 19 is T -T A .. . ,,, . , . 

& w sat As the bulk liquid is 

subcooled, the wall temperature is only slightly greater than the sat- 
uration temperature, and "AT" is therefore small. ) 

Region II is the zone of "nucleate boiling" heat transfer, and is 
characterized by the formation of bubbles at discrete points distributed 
across the heat transfer surface: i. e. , at surface imperfections, such 
as scratches, pits, etc. As the bubbles grow and break away from the 
surface, they are replaced by relatively cool liquid swept in from the 
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FIGURE 19 

"Typical" Saturated Boiling Curves for Cryogenic Liquids 
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bulk of the fluid. This bubble agitation results in a rapid increase in 
heat flux with the number of "active" nucleation sites, as the primary 
flow of heat is across the surface into the relatively cool liquid. Since 
surface imperfections tend to have a random distribution in size>, the 
number of active sites is a function of AT, with "large" sites becoming 
active before "small" sites. 

Region III represents the gradual transition from "nucleate" to "film" 
boiling. In this region, generally called the "metastable film boiling" 
region, increasingly greater portions of the surface are blanketed by 
solid "sheets" of vapor. Under such "sheets", the cool liquid cannot 
come in contact with the surface. Since the thermal conductivity of the 
vapor is generally much lower than that of the liquid phase, heat flux 
actually decreases with increasing AT, in this region. 

Region IV is the zone of fully developed "film boiling". Here, the 
entire heat transfer surface is covered by a stable vapor film, and as 
AT increases, heat flux again starts to rise. Finally, at very high AT, 
radiation effects become significant, and heat flux can become quite large. 

The heat flux attained at the transition between regions II and III 
(point "a" in Fig. 19) is variously termed the "peak", "critical", or 
"burnout" flux. The last term is quite appropriate. If heat is generated 
at a constantly increasing rate in, say, an electrically heated device, 
upon reaching this point on the nucleate boiling curve, the operating 
point must "jump" to the corresponding flux on the film boiling curve 
(point "b" in Fig. 19), since the characteristics of the device would not 
permit the decrease in flux required for operation in region III. However, 
as indicated by the "typical" scale along the abscissa, the AT required 
to maintain the same flux in region IV is generally so large that the device 
might well melt before the operating point could reach the film boiling 
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curve. 



C. Pool boiling data : 

To turn from the general to the specific, four particular cryogenic 
coolants will be considered in this study: nitrogen, neon, parahydrogen, 
and helium. All but neon have been widely used by various experimenters, 
and if current trends continue (e. g. , see (109) or (110) ), liquid neon is 
likely to be widely used in the future. As coolants, these liquids effect- 
ively span the 0°-100°K temperature range of interest here. (One-atmos- 
phere boiling points are 77. 4° K, 27. 2° K, 20. 3° K, and 4. 2° K, respec- 
tively. Of course, some control of boiling temperature is possible in 
each case by variation of operating pressures. ) 

For basic information regarding boiling heat transfer to liquid nit- 
rogen or hydrogen, one is advised to turn to either of two excellent, 
recently published works: 

Boiling Heat Transfer for Cryogenics , by J. D. Seader, W. S. Miller, 
and L. A. Kalvinskas (111), is a complete, 171-page survey of the liter- 
ature on heat transfer to liquid hydrogen, oxygen, and nitrogen. This 
work discusses and evaluates both predictive correlations and experimental 
data, for pool boiling, forced-convection subcooled boiling, forced-con- 
vection saturated boiling and forced-convection film boiling to these liq- 
uids, and includes a twenty-page bibliography, as well as eighteen pages 
of physical property charts. The work is particularly valuable for its 
coverage of the largely neglected topic of cryogenic forced convection. 

Nucleate and Film Pool Boiling Design Correlations for O 2 , N 2 , H 2 , 
and He, by E. G. Brentari and H. V. Smith (112), presents an even more 



detailed survey of pool boiling data for nitrogen and hydrogen, and, as is 
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apparent, covers helium as well. The authors discuss the various fact- 
ors which contribute most directly to "scatter" in observed data, and 
conclude that the general nucleate heat transfer correlation and maximum 
nucleate heat transfer correlation of Kutateladze (113) provide the best 
"a priori" design standards, when details of surface condition, orientation 
geometry, etc. are not known. Specifically, these correlations state: 
(q/A) = 8. 54 x 10" 1 °(c p /h lv P y ) 1, 5 °(k p 1, 282 p 1 * 750 /f 90 ^* 626 )(AT) 2, 50 
and 

<« /A) bur„out = • 157 h lvPv°' 5 °(g^>°- M <26> 
where all properties pertain to the liquid phase, except those with sub- 

i 2 

script V\ and q/A is heat flux in watts/cm ; c , specific heat at con- 
stant pressure, joules /gm-°K; h^, latent heat of vaporization, joules/ 
gm; p, density, gm/cm^; k, thermal conductivity, watts/cm-°K; P, 
pressure, dynes /cm ; t, surface tension between liquid and vapor, dynes/ 
cm; p, viscosity, gm/cm-sec; AT in °K, and g is the acceleration due 
to gravity in cm/ sec . (The bibliography provided in their work includes 
a section listing possible sources of information regarding thermodynamic 
properties. ) More conveniently, Brentari and Smith provide a comparative 
plot of these correlations and available data, and a second plot of pre- 
dicted q/A vs. AT, for various pressures (in atmospheres), for each 
liquid. 

Unfortunately, very little experimental information is available at 
this time regarding the heat transfer characteristics of liquid neon, prob- 
ably due to the hitherto limited supply. Hopefully, some data may soon 
be available as a by-product of the testing of the NASA-Lewis magnet 
(21, 22), but the only experimental report known to the author at this time 
is that of Bewilogua and Mahn (114). However, one can estimate the heat 

transfer characteristics from eqns. (25) and (26), and in fact, Frederking 
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(25) 



(ir5) does present a plot of q/A^ urnou ^, based on eqn. (26). (The National 
Bureau of Standards appears to be currently updating its information re- 
garding the thermodynamic properties of neon. See, for example, (116) 
and (117). ) 

In any event, based on (112) and (115), one can construct a quick 
estimate of the nucleate boiling heat transfer characteristics of the four 

specified coolants, as follows: Assume the one-atmosphere "burnout" 

2 o 

flux for nitrogen is 19 watts/cm , and that is occurs at a AT of 12 K. 

Similarly, for neon the figures would be 12 watts/ cm^at 3° K; for hydrogen, 
8. 5 at 3° K, and for helium, 0. 77 at 0. 4° K. Plot these points on a log- 
log chart, and pass a line downward from each, at a slope of 2. 5, to 
obtain the one-atmosphere characteristic. (Of course, since thermo- 
dynamic properties are a function of pressure, both eqns. (25) and (26) 
are functions of pressure. In general, these characteristic lines shift 
toward lower AT's as pressure is increased. The locus of the burnout 
point, as a function of pressure, is (roughly) a concave-downward arc. 

Helium exhibits its maximum burnout flux at about one atmosphere. For 

2 

the other liquids, the maximum is roughly: nitrogen, 37 watts /cm at 
AT = 6. 6°K (p = 10 atm. ); neon, 18. 6 at 1. 3°K (-'■'8. 5 atm. ), and hydrogen, 
13. 2 at 1. 6°K (4 atm. ). 

Again, it must be stressed that this is a rather crude "first estimate". 
Equation (25) is a reasonable, slightly conservative "mean", and is prob- 
ably an excellent guess for "average" conditions. However, conditions 
may deviate considerably from "average", as evidenced by the fact that 
results reported by some experimenters lie more than a factor of 5 above 
or below this correlation. For this reason, it may be worthwhile to 

briefly discuss a few of the factors which lead to deviation from the norm:' 
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Surface condition : As mentioned above, heat flux varies with the num- 

1 

ber of active nucleation sites, which among other things is a function of 
the size distribution of surface imperfections. Thus, a mirror finish 
would hinder nucleate boiling heat transfer, and data points taken from 
specimens with such a finish fall well below the norm. Conversely, an 
artificially roughened surface improves nucleate boiling heat transfer. 

(See, for example, (118). ) Further, the number of active sites is a 
function of the extent to which a particular liquid "wets" the surface. Any 
coating or contaminant which decreases "wettability" of the surface will 
tend to increase heat flux. Cowley, et al. (119) report that artificially 
greased surfaces can be driven to 50^°- 100^° higher flux levels before 
burnout, though such surfaces exhibit "below normal" flux levels at sub- 
critical AT's. Here, the dominant factor is the poor thermal conductivity 
of grease. The presence of such an "insulating layer" (~ 0. 2 mm thick) 
causes a four- to ten-fold increase in total AT (between the bulk liquid 
and the surface) at burnout. The greatest portion of this total AT occurs 
in the insulating layer, and oddly enough, results in a temperature pro- 
file which drives a greater heat flux through the insulation, into the liquid, 
than could be driven directly from the surface into the liquid. The authors 
state that such use of insulation should increase maximum heat flux in 
any situation in which heat transfer is initially limited by conditions on 
the liquid side of the interface. 

Surface history : A given surface will generally exhibit greater heat 

flux when first immersed in coolant, than after repeated measurements 

under the same conditions. This is partially attributable to the boil-off 

of surface contaminants, but is probably more dependent upon the initial 

presence of adsorbed gases. Haselden and Peters (120), having taken ( 

pains to ensure an initially clean surface, found that such transients per- 
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sisted for as long as 2 to 3 hours. Experimenters who do not choose to 
wait for such an interval, to ensure that consitions have stabilized, may 
find that their data predict higher heat fluxes than can be maintained for 
steady- state operation. 

Surface orientation and geometry : Although nucleate pool boiling is 
not particularly sensitive to device geometry, any features which alter 
the circulation pattern associated with "normal" natural convection (and 
therefore cause deviation from the strict definition of "pool boiling") will 
have some influence on nucleate boiling heat transfer. For example, con- 
fined channels may tend to become "clogged" with vapor, to such an extent 
that the churning action associated with nucleate boiling is inhibited. This 
would sharply reduce attainable heat flux. Similarly, heat transfer from 
a "bottom" surface may be hampered in a gravity field, since bouyancy 
may impair the normal motion of bubbles away from the surface. Vari- 
ations in convective circulation can also cause minor differences between 
measurements taken near the ends of a test specimen, and those taken 
near the center. However, by far the most critical of these geometric 
factors is the marked decrease in attainable heat flux in confined channels. 
Since this is precisely the type of geometry anticipated in the cooling sys- 
tems of large coils, natural convection cooling is not an attractive solution 
for such applications, 

D. Qualitative influence of forced convection : 

The obvious solution to the problem of "confined geometry" coolant 

ducts is to turn to forced convection. In addition to the "microscopic" 

processes discussed above, this adds a new form of turbulent mixing to 

the heat transfer process: i. e. , the "macroscopic" circulation of forced 

flow. Qualitatively, the effect on the nucleate boiling curve is as depicted 
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in Fig. 19. In general, the increased heat flux attributable to forced flow 

can always be made to offset or more than offset the decrement in 

attainable flux due to confined geometry. At low qualities, forced con- 
vection produces a sharp increase in both heat transfer coefficient, and 
burnout flux. As quality increases the heat transfer coefficient increases 
even more (as shown in Fig. 20a), until transition to film boiling occurs, 
after which the heat transfer coefficient decreases sharply. The forced 
convection burnout flux decreases more or less linearly with quality, as 
shown in Fig. 20b. However, even at high qualities, one might reason- 
ably expect the forced convection burnout flux to exceed the pool boiling 
burnout flux. 

The heat flux attained under forced flow conditions is often concep- 
tually divided into its two component parts: one attributable to pool 
boiling processes, and the other attributable to forced convection. Actu- 
ally, the physics of the problem are not quite that simple. As the Reynold' 
number of the flow is increased, nucleation is gradually suppressed. As 
the quality (mass fraction of vapor) of the flow is increased, the flow of 
relatively large volumes of vapor along the tube increases the convective 
turbulence attributed to forced flow. Eventually, a point is reached where 
the contribution of "microscopic" heat flux is entirely negligible, but by 
this point, the "macroscopic" heat flux may be much greater than that pre- 
dicted from eqn. (24). Seader, et al. (Ill) present a detailed discussion 
of the predictive correlations which attempt to cope with these interrelated 
complications, but in general, such correlations have not met with marked 
success, when applied to "conventional" coolants. Data for forced con- 
vection saturated boiling heat transfer to cryogenic coolants is so scarce 
that Dne can hardly say the various correlations have been fairly tested, 
for cryogenic coolants. 
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FIGURE 20a 



Qualitative Variation of Heat Transfer Coefficient 

with Quality 




FIGURE 20b 

Qualitative Variation of Critical Flux with Quality 



(q/A) 
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In any event, it would appear that the most successful such correlation, 

to date, is that developed by Chen (121): 

q/A = (h FC F + h pB S) AT, (27) 

where is the forced convection heat transfer coefficient, which could, 

for example, be computed from eqn. (24); hpg is the pool boiling heat 

transfer coefficient, which could be computed by dividing eqn. (25) by AT; 

F is an empirical amplification factor which corrects for increased heat 

flux due to increased convective turbulence at high quality, and S is an 

empirical attenuation factor, which corrects for the suppression of 

nucleation at high liquid flow rates. Both (111) and (121) present plots 

of F as a function of the Martinelli two-phase parameter, X^, and S as a 
1 ?5 

function of R T F ' , where 

eL ’ 



lO. 9 



tt 



\A ] 

Pi 



\0. 5 




0. 9 



(28) 



(here, x is quality, and the other parameters are as defined above), and 

R eL = DG<1 - x)/ Ml <29) 

(where G is the mass flow rate divided by the tube cross sectional area). 
As a practical matter, the plots provided in (111) are, perhaps, somewhat 
easier to work with. 

It is apparent that the use of this correlation to predict heat flux at a 



single point along a coolant duct is no minor task, unless the problem has 
been prepared in advance for computer solution. Further, since both 
eqn. (28) and eqn. (29) are functions of quality, it is necessary to evaluate 
eqn. (27) at many points along such a duct, if a "mean" flux is to be deter- 
mined for purposes of preliminary analysis. Chen has tested this cor- 
relation against available data for vertical, axial, stable flow, at heat 
fluxes below the critical, for water and various hydrocarbons, and found 
an average deviation of only ±11 percent, for a quality range of 1 to 71 
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percent. This achievement is truly remarkable, considering the com- 
plexity of the problem! However, individual points may deviate consid- 
erably more than the "average", and forced convection heat transfer is 
more sensitive to geometry and orientation than pool boiling heat transfer. 
Considering this, and the scarcity of data for two-phase forced convection 
heat transfer to cryogenic coolants, it should be stressed again that experi- 
mental confirmation of heat transfer characteristics should be obtained 
before proceeding into the detailed design phase. 

Under the above circumstances, it is doubtful that the considerable 
labor involved in the use of Chen's correlation is justifiable for purposes 
of preliminary analysis. Instead, it is recommended here that the pro- 
cedure for determining pool boiling heat transfer characteristics, as 
described in the preceding section, be applied with no correction for 
forced flow. This procedure has the advantage of being quicker, simpler, 
and conservative, since the forced convection heat transfer characteristics 
should be as good as, or better than, this estimate. (It is noted that 
Frederking (122) specifies that heat flux under conditions of forced con- 
vection will always exceed predictions from a pool boiling analysis if 
quality is low . However, Chen's correlation is less restrictive than this, 
since fluxes predicted from his model exceed the pool boiling figures up 
to the limiting quality of 70%. ) 
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VI. COOLANT FLOW 



A. Ge ometric considerations: 

While heat flux considerations may determine the amount of heat trans- 
fer surface required for a particular coil, the specification of cooling sys- 
tem geometry is far from complete. The required surface area could be 
provided in a single large coolant duct, or many small ducts, and the choice 
between these extremes must be based on considerations other than heat 
flow alone, such as the resulting packing factor, pressure drop across the 
coil, thermal gradients in the winding, ease of assembly, etc. , etc. 

To deal with the last item, in particular: The natural tendency would 
be to assume some specified distribution of circular tubes, of some given 
L/D. This would be motivated in part by historical precedent, and in part 
by convenience: Most conventional heat transfer devices utilize circular 
tubes; we are accustomed to thinking in terms of circular tubes, and the 
great majority of the correlations used in heat transfer and coolant flow 
analyses are based on circular tube data. Unfortunately, circular tubes 
are not likely to be used in "real" cooling systems for coils of the type 
being considered here, for precisely the same reason that they are used 
in so many other applications: i. e. , ease of manufacture. 

Whatever provisions are made for coolant flow must conform to the 

realities of coil construction. For example, it was shown in Chapter III 

that the internal support structure required to support the conductors 

against magnetic stress must run axially down the length of the cylindrical 

lobes. Thus, it is reasonable to assume that the coolant ducts must run 

axially: any other geometry would require the ducts to penetrate such 

structural members, which would be highly undesirable, if not impossible. 

Further, cooling ducts cannot be run through regions in which they would 
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be subjected to a heavy stress load. This would require that they have 
thick walls, for mechanical strength, which would extract a double pen- 
alty: reduced heat transfer, and reduced packing factor. Since structure 
has already been provided to support the conductors, it is reasonable to 
assume that coolant ducts would be located contiguous to that structure: 
stress levels are at a maximum on one side of each support shell, but 
total stress is virtually zero on the other side. A coolant duct, run along 
the "no-load" side of such a shell, would require a minimum of structural 
support. This fact is conducive to the use of coolant slots , vice tubes, 
or pipes. Of course, the model described in Chapter III calls for a greater 
number of support members in regions of high magnetic fields, since these 
are the regions of greatest magnetic force density. These are also the 
regions where magnetoresistance, and therefore heat generation, would 
be greatest. It is convenient that a greater number of coolant slots would 
be provided automatically in this region, under the proposed scheme. 

Thus, the basic geometry considered here will involve the use of 
axial slots for coolant flow, with one side of each slot bounded by the 
internal support structure of the coil. It will be assumed that the height 
of such slots would be much less than all other significant dimensions 
(width, length, and radius of curvature), and coolant flow will therefore 
be handled on a "per unit width" basis. Since the most likely structural 
material is stainless steel, and since the thermal conductivity of stain- 
less steel is two to four orders of magnitude less than that of the electrical 
conductors considered here, at the temperatures of interest, it will be 
assumed that heat enters each coolant slot only across the surface "away 
from" the support shell. (So far as it involves the results of this analysis, 
one can treat the possibility of heat input through both surfaces by simply 

doubling the attainable heat flux which would otherwise be estimated from 
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the preceding chapter. ) Again, the purpose of this analysis will be to deter- 
mine the relationships between coolant properties, pressure, quality, and 
heat flux, as "input" variables, and packing factor, coolant pump power, 
and axial temperature variation, as "output" variables. 

B. P roperties of coolants : 

Since coolant properties are a necessary input to this analysis, it 
would be well to consider them before proceeding further. In particular, 
since the primary effort here will be devoted to an analysis of two-phase 
cooling, the saturation properties (i. e. , the thermodynamic properties of 
a saturated liquid and its vapor, at least nominally under equilibrium con- 
ditions) are of especial interest. In fact, an interesting problem exists 
regarding the processing of this data. Because the coolants considered 
(helium, parahydrogen, neon and nitrogen) will be examined at pressures 
approaching their respective critical pressures, such customary and con- 
venient concepts as "incompressible" liquids, and "perfect" gases, are 
not generally applicable. Further, many pertinent properties exhibit 
sufficiently non-linear behavior, within the range of interest, as to raise 
doubts regarding the accuracy of such techniques as linear interpolation 
between tabulated values. As a convenient means to avoid these compli- 
cations, a data fitting routine was used to generate a least -square poly- 
nomial fit, expressing each such property of interest as an explicit func - 
tion of pressure, in atmospheres. Specifically, the properties of specific 
volume of saturated liquid and vapor, saturation temperature, specific 
enthalpy of the saturated liquid, and heat of vaporization, were expressed 
as polynomials of the below form: 

2 3 

!/p i = a i + a 2P + a 3P + a 4P 

ll P v = b l /p + b 2 + b 3 p + b 4 p2 + b 5 p3 + b 6 p4 
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TABLE 6 



SATURATION PROPERTY COEFFICIENTS 



Property 


Helium 


Parahydrogen 


Neon 


Nitrogen 




f a l = 


5. 45596 


13. 3964 


0. 767861 


1. 18952 




r a 9 = 


5. 09263 


0. 796955 


0. 0675544 


0. 0559020 


fa T = 


3. 69782 


-0. 0606407 


-0. 00664262 


-0. 00437871 




a? = 


1. 15685 


0. 00500739 


0. 000224027 


0. 000195317 


Deviation = 


0. 058% 


0. 041% 


under 0. 01% 


0. 074% 




^1 = 


49. 3439 


654. 454 


101. 175 


195. 349 




b 2 = 


31. 7682 


128. 968 


0. 394318 


29. 5132 


i/Pv 


b ? = 


-32. 7486 


-41. 6982 


3. 23665 


-8. 96109 


b 4 = 


15. 3537 


6. 73209 


-0. 949753 


1. 49260 




b 5 = 


-4. 01758 


-0. 563735 


0. 0978116 


-0. 126834 




lb 6 = 
tion = 


0. 281743 


0. 0181262 


-0. 00343717 


0. 00420598 


Devia 


under 0. 01% 


0. 100% 


0. 561% 


0. 059% 




( C 1 = 


2. 67907 


16. 4517 


23. 1925 


67. 7503 




C 2 = 


2. 04654 


4. 71658 


5. 04433 


11. 9973 


T ( 


C 9 = 


-0. 461891 


-1. 00566 


-1. 25994 


-2. 66612 


o 

C A = 


-0. 180713 


0. 144354 


0. 207 938 


0. 392040 




C c = 


0. 158552 


-0. 0108406 


-0. 0176176 


-0. 0300996 




II II 

d 

n coo 
o 


-0. 0309227 


0. 000323253 


0. 000582154 


0. 000918179 


Devia 


under 0. 01% 


0. 049% 


0. 038% 


0. 035% 




'd 1 - 
do 


-28. 9212 


-293. 057 


-2. 09216 


9. 57010 




139. 468 


43. 7316 


8. 43141 


24. 7498 




d 3 = 


-211. 281 


-7. 27746 


-1. 56955 


-5. 51273 


d 4 = 


163. 231 


1. 05418 


0. 201394 


0. 820578 




d 5 = 


-61. 7927 


-0. 0791551 


-0. 0130336 


-0. 0632501 




idl = 


9. 19968 


0. 00245513 


0. 000326842 


0. 00193276 


Deviation = 


under 0. 01% 


0. 033% 


0. 721% 


0. 134% 




' e 1 = 


58. 3544 


462. 481 


91. 6010 


214. 510 




e n = 


-137. 426 


-14. 4022 


-7. 50772 


-18. 9988 


h LV< 


e 3 = 


211. 270 


-1. 53707 


1. 64718 


4. 77690 


I e 4 = 


-163. 972 


0. 159357 


-0. 200480 


-0. 789073 




e 5 = 


61. 7927 


-0. 00838393 


0. 0130336 


0. 0632501 


1 


3® = 


-9. 19968 


0. 00245513 


-0. 000326842 


-0. 00193276 


Deviation = 


0. 035% 


0. 015% 


0. 693% 


0. 148% 
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T = C 1 + c 2 p + c 3 p 2 + c 4 p 3 + c 5 p 4 + c 6 p 5 
h l = d i + d 2 p + d 3 p + d 4 p J + d 5 p + dgp^ 

h lv = e l + e 2 P + e 3 p2 + e 4 p3 + e 5 p4 + e 6 p5 
where specific volumes are given in cm 3 /gm, temperature in °K, and 

enthalpies in joules /gm. The coefficients a., b^ , etc., are tabulated 
below. The particular functional forms indicated ensure, in all cases, that 
the accuracy of "fit" exceeds the assumed accuracy of the flow models 
developed here, within the pressure ranges investigated: i. e. , 1-10 atmos- 
pheres for saturated nitrogen, neon, and parahydrogen, and 1-2 atmos- 
pheres for saturated helium. (The lower limit of 1 atmosphere was set 
because it was felt that considering the hazard of leakage, no one would 
care to design a cryogenic cooling system, for a working coil, to operate 
at sub-atmospheric pressures. The upper limits were set, as a practical 
matter, considering the decrease in heat of vaporization, and increase in 
coil operating temperature, at increased operating pressures. ) The "devi- 
ation" indicated for each set of coefficients is the mean of the absolute 
value of the difference between "input" property values, and property val- 
ues generated from the indicated polynomial. 

One other general point might be mentioned before turning to the par- 
ticular: If two-phase cooling is contemplated, one would presumably 
attempt to operate to a relatively high outlet quality, since pumping "unused" 
liquid through the cooling loop increases pump capacity, and increases duct 
size, while contributing nothing to the removal of heat from the coil. In 
practice, outlet quality must be limited by heat transfer considerations 
(more on that, later!), but if attainable qualities even approach the range 
indicated by the upper limit of Chen’s correlation, it should be recognized 
that the outlet stream, from a volumetric viewpoint, will be almost entirely 

vapor: a "typical" figure would be 99% vapor, for a quality of 50%. As 
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vapor densities tend to be relatively low, momentum considerations indi- 
cate that relatively high vapor outlet velocities may result from relatively 
low pressure drops: i. e. , it is quite possible that the flow could "choke" 
at, or even before, the outlet. If this is to be avoided, outlet velocities 
must not be allowed to approach the sonic velocity in the stream. (Note 
that this would be the so-called "sonic velocity of Mach l", which Shapiro 
(123) denotes "c*", not the velocity of sound under stagnation conditions. ) 
The problem here is to determine just what the limiting velocity should be. 
The only available data are for c , the sonic velocity under stagnation 
conditions. Shapiro shows that for perfect gases: 

c*/c Q = V2/(7+ 1), 

where = c /c = the ratio of specific heats. However, we are not deal- 
p v 

ing with perfect gases, and in any event, the sonic velocity in the mixed- 
phase coolant stream would presumably be lower than that in the vapor, 
under conditions of 100% quality. Considering these uncertainties, a 
"limiting velocity" has been defined for each coolant, by using the best 
available information for c q and ^ at the one atmosphere saturation point, 
computing c* as above, and multiplying the result by roughly 0. 8. In the 
subsequent analyses, combinations of parameters which would result in 
flow velocities greater than this limiting velocity will not be considered. 
The values of , c q , and v^ are presented below, along with other spe- 
cific information of interest for each coolant. 

Helium : (Property values determined form (102,110) and (124-127). 

Saturation property polynomials developed from (127). ) 

one -atmosphere boiling point: 4. 22°K 

critical point: 5. 25 K at 2. 26 atm. 

ratio of specific heats: ~1. 74 

sonic velocity: ~145 m/s 

limiting velocity: 100 m/s 

Considering its low heat of vaporization, and the high cost of 4. 2°K 
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refrigeration, liquid helium is not likely to be used as a two-phase coolant 
for cryogenic coils, except when cryogenic conductors are used to "stabi- 
lize" superconducting coils. (See, e. g. , (14). ) However, supercritical 
helium is the logical coolant for coils designed to operate in the vicinity 
of 10°K (24, 107). The National Bureau of Standards indicates that the 
optimum operating pressures for such applications lie in the band 15-20 
atm. (105). The required coolant volume fraction is not a strong function 
of pressure, in this range, since the increase in density as pressure is 
increased from 15 to 20 atm. is largely offset by a decrease in specific 
heat. 

Parahydrogen: (Property values determined from (102, 110, 124-126) 

and (128-131). Saturation property polynomials developed from (131), 

except that for liquid specific enthalpy, which was developed from (130). ) 

one-atmosphere boiling point: 20. 27°K 

critical point: 32. 98°K at 12. 76 atm. 

ratio of specific heats: 1. 803 

sonic velocity: 357 m/s 

limiting velocity: 250 m/s 

Hydrogen gas occurs in two molecular forms, orthohydrogen and para- 
hydrogen, distinguished by the relative direction of nuclear spin of the two 
component atoms of the molecule. At room temperature, the equilibrium 
composition is 75 per cent ortho, 25 per cent para. However, at 20. 3°K, 
the equilibrium concentration is about 99% parahydrogen. Since the energy 
released during conversion is greater than the latent heat of vaporization, 
this reaction must be catalyzed to completion if liquid hydrogen is to be 
used continuously in closed-loop cooling systems. Presumably, the com- 
position would not change significantly during the brief period required for 
the liquid to pass through the coils of a magnet. For these reasons, the 
properties discussed here are those of liquid parahydrogen. 
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Basically, liquid hydrogen is a highly promising coolant for cryogenic 
coils, because of its availability and low boiling point. Although some con- 
cern has been expressed regarding the safety of this particular coolant, 
those who have worked with liquid hydrogen (notably, (19, 20) report "no 
problems", provided reasonable precautions are observed. (See NBS 
Report 7253 (132) for a list of safety precautions. ) It would seem reason- 
able to presume that the hazards involved in the use of this coolant in 
closed-loop systems would be less severe than is the case with open-loop 
systems, where large volumes of hydrogen gas must be stored or vented. 

Neon : (Property values determined from (102, 110, 133) and (116). 

Saturation property polynomials developed from (116). ) 

one-atmosphere boiling point: 27. 09°K 

critical point: 44. 40°K at 2 6. 19 atm. 

ratio of specific heats: ~1. 7 

sonic velocity: 134 m/s 

limiting velocity: 90 m/s 

As has been discussed above, liquid neon has not been widely used as 
a cryogenic coolant, to date, because of its hitherto limited availability 
and high cost. Specific advantages attributed to this liquid as a coolant 
are a relatively high latent heat, on a per-unit-volume basis, and ease of 
temperature control through pressure regulation. In addition, liquid neon 
is the only two-phase coolant available which permits maintenance of tem- 
peratures nearlthe liquid hydrogen range, without the hazards associated 
with liquid hydrogen. 

Nitrogen: (Property values determined from (102, 124-126, 128, 129) 

and (134). Saturation property polynomials developed from (134). ) 

one -atmosphere boiling point: 77. 4°K 

critical point: 126. 1 K 

ratio of specific heats: 1. 44 

sonic velocity: ~167 m/s 

limiting velocity: 120 m/s 

Nitrogen might be described as the "old standby" cryogenic coolant, 
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and apparently is used in many applications primarily because its proper- 
ties are well known, it is readily available in large quantities, and it can 
be handled with relative ease. However, although various "laboratory" 
coils have been tested in liquid nitrogen. Tables 1-4 indicate that the tem- 
peratures attainable with this coolant are too high to permit realization of 
the full potential of cryogenic magnets. It is unlikely that large "working" 
coils would be designed for two-phase liquid nitrogen cooling. 

C. Single-phase flow : 

Aside from the concepts discussed in the preceding chapter, the only 
design factor essential to the preliminary analysis of single-phase cooling 
is the question of pressure drop across the coil. The customary manner 
of estimating this pressure drop would be to do so in terms of an effective 
friction factor, f: 

p = 4 f (L/D)(pv 2 /2), (30) 

where one would again utilize "mean" property values, and where L is the 
estimated duct length, which would presumably be of the order of the MHD 
channel length (or, perhaps, the product of that length times the form fac- 
tor, defined in eqn. (6) ). For non- circular cross sections, D would be 
the "equivalent" or "hydraulic" diameter, given by: 

D = 4 A/P (31) 

where A is the duct cross-sectional area, and P is the wetted perimeter. 
(For a rectangular slot, of the type discussed previously, 4A-*-4ws and P 
-*-2w, so D S' 2s, since slot height s » slot width w. ) 

Various estimates are available for the value of f, given fully developed 

I 

turbulent flow conditions. Rohsenow and Choi (108) suggest: 

f = 0. 0791/R e 0, 25 , 3 x 10 3 < R 0 =2 x 10 4 (32) 
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or 



(32) 



f = 0. 046/R e °' 2 °, R g >2x 10 4 . 

Stevenson (107), on the other hand, uses:' 

0 19 

f = 0. 0014 + 0. 125/R * . (33) 

(In both eqns. (32) and eqn. (33), Rg-pvD/^/.) 

Having established an estimated bulk flow rate (from mean temperature 

rise considerations), the above Ap can be used to determine the power 

required to pump the coolant through the loop, since 

P ? Fv ? (AAp)v - wAp/p (34) 

where w is now the mass flow rate. Of course, this term represents 

energy added to the coolant stream, which is converted to heat by friction, 

and is therefore an added refrigeration load. However, one would antic- 

2 

ipate that this load would normally be negligible, compared to I R heating. 
Eqn. (34) can be used to verify this fact, as well as to "size" the pump 
required for a particular coil. Once again, the procedure is relatively 
straightforward, and a great deal simpler than the analysis of two-phase 
flow. 

D. Two-phase flow : 

Considering the lack of experimental data for forced- convection two- 
phase cryogenic cooling, the information required for coolant flow calcu- 
lations must be generated from an analytical model. The problem becomes 
to select that model which is most likely to correspond to the physical 
realities of the situation, under these circumstances. Two possibilities 
will be developed and discussed here, but the ultimate choice must await 
additional testing. 

To turn, first, to background considerations: Sydoriak and Roberts 

(128) conducted one of the few studies of two-phase cryogenic coolant flow 
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in confined channels, for both natural and forced convection, and report 
close correlation between their data and "calculations based on a friction- 
less flow boiling theory". In essence, the pressure drop due to the always - 
low viscosity of cryogenic fluids was totally dominated by that attributable 
to acceleration of the vapor phase (as required to maintain a given flow 
rate when a quantity of relatively dense liquid is converted to relatively 
light vapor). The pressure drop associated with any such "inviscid" accel- 
eration can be computed from momentum considerations, provided certain 
basic information is available regarding flow conditions in the duct. 

Chen (121) classifies the possible flow conditions, as depicted in Fig. 

21 : 

Zone 1, at the inlet, represents subcooled conditions, as would be 
encountered when the saturated liquid, from a heat exchanger, passes 
through a pump, with resultant increase in pressure, but little increase 
in temperature. In this region, there is either no generation of vapor, or 
formation of bubbles, at the duct wall, which decrease in size as they depart 
the thermal boundary layer. 

Zone 2 represents the region of first net generation of vapor. Here, 
bubbles are formed at the heat transfer surface, break away, and move 
into the (now saturated) liquid phase. Quality is low, and aside from the 
convection down the duct, conditions are rather similar to "nucleate pool 
boiling". 

By zone 3, sufficient vapor has been generated that the individual 
bubbles tend to merge, forming relatively large "blobs" of vapor which 
accelerate down the duct. Nucleation is gradually suppressed. This flow 
condition is characterized as "slug flow". 

Eventually, the center of the duct becomes filled with an unbroken 

stream of vapor, as in zone 4, and an effective separation of the phases 
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'Typical" Flow Development 



FIGURE 21 
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occurs. As Chen's work dealt with circular tubes, he called this the 
"annular flow" regime. Although some liquid is carried with the vapor 
in mist-like droplets, the bulk of the liquid phase moves along the walls 
of the tube. 

Finally, as quality increases, the tube walls gradually become dry 
and zone 5, the "liquid deficient zone", is entered. All the remaining 
liquid, in this zone, is swept along in the form of droplets, dispersed 
through the vapor phase. 

Actually, the scale of Fig. 21 is distorted. Chen's correlation was 
based upon the assumption of annular flow conditions, and as mentioned 
above, is held to be valid, generally, within the range 1 to 71 per cent 
quality, for the vertical, circular-tube geometry which he investigated. 
However, the fact that he determined this to be true for one geometry is 
no assurance that it will be true for other geometries . The various re- 
gimes may be less clear-cut, with some zones missing, others enlarged, 
or mingled. The sequence depicted here is "idealized", and deviations 
from this sequence are important in two respects: first, flow conditions 
in zone 4 are most nearly typified by separate, inviscid flow of the liquid 
and vapor phases, with no mixing between the phases. Zone 5, and per- 
haps zones 1 and 3, are better typified as inviscid flow, with the boundry 
condition of "no slip" between the phases. Thus, any shift in the relative 
importance of these zones would presumably require a shift in the analyt- 
ical model used to represent flow conditions. Second, when the tube walls 
become dry, the heat transfer coefficient drops abruptly, due to the low 
thermal conductivity of the vapor phase, which then contacts the heat trans- 
fer surface. If zone 5 is entered at qualities significantly below 70% , for 
a particular geometry, the outlet quality for a coil cooled by such a system 
would perforce be limited to such lower qualities. In other words, the fact 
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that Chen has proven that heat transfer improves with quality if annular 
flow conditions are maintained, does not, a priori, guarantee that such 
conditions will be maintained. Again, conditions must be confirmed for 
the particular application contemplated. 

Bearing this in mind, two inviscid flow models will be examined 
below: 

No slip (mist flow) model : Of the two limiting conditions discussed above, 
the simpler, from an analytical viewpoint, is that assuming no slip between 
the phases. If either phase is assumed to be dispersed through, and moving 
with, the other, the resultant two-phase flow can be treated as a hypo- 
thetical, homogeneous fluid with properties which are appropriately aver- 
aged between those of the liquid and the vapor. Specifically, let us assume 
steady, inviscid, one-dimensional flow, traversing a slot of constant, but 
as yet unspecified height, of unit width, and of some length L. The vari- 
able z will be used to denote distance from the inlet of the slot (x will refer 
to quality). It will be assumed that the flow has just reached zero quality 
at the slot inlet (i. e. , x=0 at z = 0), and is subjected to a constant heat flux, 
q/A, between z = 0 and z = L. (This last is equivalent to the assumption that 
heat transfer to "real" liquids, which is known to be a function of quality, 
can be expressed in terms of an effective "mean" q/A. ) With an eye toward 
eventual applications, it will be assumed that outlet pressure, and desired 
outlet quality, are specified, and that we desire to investigate the relation- 
ship between inlet pressure (or, equivalently, pressure drop across the 
slot), inlet velocity, and slot height. The axial distribution of temperature, 
velocity, etc. will also be of interest, and in particular, the solution for 
each combination of parameters must be tested to ensure that the outlet 
velocity does not exceed the limiting velocity, defined in Section B of this 
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An iterative approach will be used to accomplish these objectives. In 
order to investigate the influences of inlet pressure, that parameter will 
be systematically incremented, as: 

P t = P e + n Ap n = 1, 2, 3 

By use of the saturation property polynomials, all pertinent thermo- 
dynamic properties at the inlet and outlet can be determined, once p^ and 
P Q have been specified. However, the "appropriately averaged" density 
of the "homogeneous" fluid is required, to relate these variables to flow 
rate. To quantize this concept: one gram of coolant, in the duct, is by 
definition comprised of x grams of vapor, and (1-x) grams of liquid. The 
vapor occupies a volume of x/ cubic centimeters, and the liquid, (1-x) 
/ P y Since the total volume occupied by this gram of coolant is, by defi- 
nition, the specific volume of the hypothetical fluid, the "appropriately 
averaged" density is defined as: 

P - */p v +»-x)/ Pl < 35 > 

Since we are considering a slot of unit width, continuity implies: 
w = p-v A = pvs = const. , (36) 

where w is the mass flow rate, v the common velocity of the two phases 
(or of the "homogeneous" fluid), and s the still undetermined slot height. 
However, since s is constant, 

pv = const. = jQ.v^ (37) 

(here, as elsewhere, subscript "i" refers to inlet properties). Under the 
stipulated flow conditions, the momentum equation can be integrated di- 
rectly, between inlet and outlet, to yield : 

-Ap = p Q - = ( P i v i ) Av = (PiV^/Po " 1/ Pi ) (38) 

Since the only unknown in the above is v^, we can immediately state: 
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Vj = (l/p i )VIArf/(l/p 0 - 1/p t ) 



(39) 



The same relationship can also be solved for v q , which can immedi- 
ately be checked against the specified limiting velocity: 

v o = v i+ P 1 i v i (40) 

Assuming v q < v^ , we proceed: Since inlet quality was stipulated 
to be zero, and the slot is of unit width, at a point infinitesimally removed 
from the inlet, by a distance dz, we have: 



_ (q/A) dz 



w h 



(41) 



lv 



Here, we encounter two problems. First, although h^ is a known function 

of pressure, the pressure distribution along the slot has not yet been 

determined, so eqn. (41) cannot be integrated directly to determine x as 

a function of z. Second, though the ratio (q/A)/w is known to be constant, 

s is not known in eqn. (36), and the value of this constant is therefore 

unknown. Letting the unknown constant be K, one additional item is^ known: 
L 



•'rv 



(q/A) dz 



w h 



lv 



where 



lv, 




(42) 



Of course, (1/h^) is not known, either, but the problem is now in 
convenient form for an iterative, numerical solution. Assume 



lv 



1 _ 

2 



h, h, 
lv i lv o< 



(43) 



By use of eqn. (43), eqn. (42) may be solved for a "first guess" for K. 

Divide the slot into N segments, of length Az = L/N. If N is sufficiently 

large (N ^10 yields reasonable results), one may assume^ 3 Pi ,h lv = blv^ 

etc. , throughout the first segment. Eqn. (41) may then be used to determine 

the quality at the end of the first segment, after which eqn. (35) can be used 

to increment p , then eqn. (37) to increment v, and eqn. (38) to increment p. 

Ill 



With p incremented, new values may be determined for the thermodynamic 
properties, and these used to advance the solution through the next segment, 
etc. When the outlet is reached, computed outlet quality can be compared 
to desired outlet quality, and the estimate for K corrected accordingly. 

A simple, linear correction was found to yield rapid convergence: 

K = K.(x, ,/x ,), etc. 

2 1 desd compt ' 

Of course, when suitable convergence is obtained, the desired tem- 
perature, pressure, and velocity profiles have also been obtained, and 
eqn. (42) may be solved directly for the hitherto unknown relationship 
between (q/A) and s. 

No mixing (annular flow) model : Essentially, the assumptions made regard- 
ing flow conditions for the no- slip model pertain as well to the no-mixing 
case, except that concerning the distribution of the phases. Here, we 
shall assume that the liquid phase remains adjacent to the walls of the slot, 
while the vapor phase flows along the center of the slot. With neither phase 
entrained in the other, the assumption of inviscid flow implies that liquid 
and vapor can (and indeed, must ! ) proceed at different velocities. In- 
stead, if the pressure profile is to be one-dimensional in character, the 
bulk velocity of each phase must adjust (considering density differences) 
so that the pressure drop along, say, a streamline in the center of the duct 

which lies entirely in the vapor phase is equal, at any z, to the 

pressure drop along a streamline adjacent to the wall, lying entirely in 
the liquid phase. Thus, we now have two separate momentum equations to 
deal with, and therefore two different phase velocities, which precludes 
treating the two-phase flow as a simple "homogeneous" fluid. (At the inlet, 
assuming quality is exactly zero, the first infinitesimal quantity of vapor 

generated causes no finite change in flow conditions. Therefore, at the 
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inlet, and only at the inlet, the two phase velocities are equal. ) 

The concept of "quality" becomes somewhat more complex, given 
that liquid and vapor velocities are generally unequal. Actually, quality 
is a static concept, which was easily extended to the preceding dynamic 
case, because the phases moved together. To avoid confusion, the term 



will here be taken to mean: 



.Ox. 



s dz 



£v S v 



O s dz + As, dz D s + O. s. 

r-y v n l ~v v r'l 1 



(44) 



where we again speak in terms of unit width, and s v is that portion of the 

slot height occupied by the vapor phase (s v + s^ = s, as defined above). 

For convenience, the variable £ will be defined as: 

. A, Si 1-x D 

t'lr-T'—TT <«> 

V V ~1 

In terms of £ , the continuity equation, written between inlet and outlet, 
now becomes: 

A. v i ( i + £ > " Pi v i £ + Pv v v (46) 

1 o o° o o 

Because specific volume is a known function of pressure, the two 

momentum equations can be integrated directly between inlet and outlet: 

Ap, 



A p v =i 



where 



P/Py) J 

Tj = fa'P 



. 2 2 X 

(v i - v i > 
o 


(47) 


, 2 2i 

(v - V. ) 

v i 


(48) 



Pe "Pi 



which can be evaluated by direct integration of the appropriate saturation 

property polynominals. Of course, as was stated above, we require 

A P 1 = Ap y (49) 

In principle, eqns. (45-49) can be reduced to three equations with the 

three unknowns v^, v , and could therefore be solved explicitly for the 
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inlet and outlet velocities. In practice, this is messy, and an iterative 

approach was employed here. This problem was solved in tandem with 

the no-slip problem, preceding. The inlet velocity determined for the 

former problem was inserted in eqns. (47) and (48), which were then 

solved for v, and v . These values were inserted in eqn. (46), which 
J-o v o 

was solved for a corrected "guess" for v., etc. Convergence was found 

to be quite rapid: generally, the value of v converged in a single iter- 

v o 

ation, and v. and v, in two or three. 
io 1 

Of course, once determined, v v can be checked against v ^ m * and 
assuming v v < we again proceed. With inlet and outlet variables 

completely determined, it is again convenient to adopt an iterative, numer- 
ical solution. In this case, the energy relationship becomes: 

(q/A) dz = h^ v (5 W V , (50) 

where (5 W V is the change in the vapor mass flow rate while traversing the 
distance dz. For an incremental analysis, let 



K 



= (q/A) L 



Ns 



(51) 



Again, s is unknown, but if the two problems are solved in tandem, we can 
use the slot height determined from the no-slip model for the same set of 
thermodynamic parameters. Using an approach similar to that of the 
preceding problem, at the k** 1 section: 
w /s=w /s+K/h. . 

v k Vi k -‘ 

Specifically, at the inlet, w v is defined to be zero. Further, w v + w^ = w, 

which is still constant. To advance the solution from one station to the 

next, let s^ = s^/s and Sj 1 = s^/s ( s v ' + s^' = 1). Then, 

v s = (1 / Q )(w /s+K/h, ) = v (1-s. ') (52) 

v k v k v k-i lv k-i v k + 



Vl' = (1/ Pi. >< w l. , /s - K/h lv. > 

k k k— 1 k-1 k-1 



(53) 
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( 54 ) 



and, to keep Ap^ = Ap y between stations: 

A, <\ 2 - v, 2 ) ■ P„ (v 



2 2 
- v 



L k-1 "k "k-l ' v l-l v k v k-l 



) 



At the inlet, s v = 0, s^ = 1, w v = ), and w^/s = Eqns. (52-54) 

can be used to advance the solution consecutively from one section to the 

next, although an iterative solution is again convenient. Here, eqn. (53) 

was solved for s.', under the assumption that 
x k 

V i " V l, . 

~ , o k- 1 

v. = v, + , 

k k-1 N-(k-l) 



eqn. (52) was then solved for v , after which eqn. (54) was used to derive 

v k 

a correction for the assumed value of v, . 

*k 

At the outlet, computed and desired outlet qualities should be compared, 
and the process repeated with a new value of K, until desired convergence 
is obtained. As was the case previously (when eqn. (42) was used to com- 
pute s), eqn. (51) cam be solved for that variable, after the solution has 
converged. 



Utilization of the above models: actually, both these problems were pro- 

grammed for computer solution, but in terms of the parameter S, which 
is defined by the relationship: 

S = s/(q/A)L (55) 

The reason for this shift of variables is that the slot factor, S, is the more 
convenient parameter for the contemplated application: for a specified p Q 
and x q , a single curve can be developed, giving S as a function of Ap. A 
family of such curves, for various x q , concisely displays the general 
relationship between inlet and outlet pressure, quality, and slot height, 
for any (q/A) and L. 

Further, while the incremental models give access to a temperature 

profile, through application of the saturation temperature polynomial to the 
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-computed pressure at each station, a more useful figure for preliminary 
analysis would be the mean coolant temperature, or better still, the mean 
temperature increment. The former is defined as: 

(5s) 

k= 1 _ 

The latter might be defined as AT = T - T , but as a matter of conven- 

c c o 

ience, will be defined instead as: 

AT = T - T , . , ,,„v 

c c sat, 1 atm. (57) 

under the assumption that the one-atmosphere saturation temperature may 

be more readily available than the outlet temperature, if outlet pressure 

differs from one atmosphere. 

The primary advantage of dealing with the temperature increment, 
defined by eqn. (57), rather than the mean temperature, is that the mag- 
nitude of the increment clearly underlines the added "cost" of operating 
at elevated pressures. In computing mean coil resistivity, the first esti- 
mate for mean coil temperature would be given by: 

T~T = T . + AT + AT,, (58) 

w s, 1 c bl 

Here, T^ is actually the wall temperature of the cooling duct, and AT^^ 
is the temperature increment across the thermal boundary layer, which 
would presumably be determined from the same source as the figure used 
for mean (q/A): either a predictive correlation, as discussed in the pre- 
ceding chapter, or experimental data. This is classed as a "first estimate", 
since if the coolant ducts are widely spaced, an additional correction might 
be necessary to account for thermal gradients in the conductor material. 
However, this factor is again best left for the detailed design phase. In 
general, the T defined in eqn. (58) would be the temperature used to enter 
the mean coil resistivity curves of Appendix B. 
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Correction for subcooling at inlet: Either of the above models could be 



applied to the problem at hand e. g. , to predict coolant volume fraction 

for a cryogenic coil depending upon which was most representative of 

l 

flow conditions realized in the type of slot geometry contemplated. Again, 
it is to be regretted that experimental data are simply not available for 
forced flow, two-phase cryogenic cooling with this type of slot geometry. 
Presumably, all efforts would be made to maintain operation in the ''annu- 
lar flow" regime, in which case the no-mixing model would seem most 
pertinent. However, in the absence of experimental data, the no-slip 
model has the advantage of producing a more conservative prediction: 

Since the relatively dense liquid must accelerate with the vapor phase in 
that case, flow velocities tend to be lower for a given Ap. This, in turn, 
requires larger slots to pass the same bulk rate of flow, as would be 
required to maintain the same specified outlet quality. Until experimental 
data confirm that viscous effects are truly negligible, throughout the range 
considered here, this "built-in conservatism" might well be used as "ad- 
vance compensation" for possible viscous effects. 

In any event, before either model can be applied to predict flow rela- 
tionships in "real" coolant ducts, one additional correction should be intro- 
duced. Both models assumed zero quality at the slot inlet, but as was 
mentioned above, the zero-quality saturation point would not, in fact, be 
reached until the subcooled liquid coolant, coming from the pump(s), had 
been heated to its saturation temperature. Particularly at high pressures, 
where h^ v tends to be relatively small, the energy absorbed before reaching 
the saturation point might be significant. (Note, however, that it should 
not be predominant, or this could not properly be classed as "two-phase 
cooling". However, the pressure ranges selected for the various coolants 
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in effect constrain the problem to avoid that possibility. ) 



What the above models actually predict is the relationship between 
pressure, slot factor, etc. , and the length of that portion of the slot 
devoted to two-phase flow. Suppose zero quality is reached at the point 
z , rather than at the inlet. The liquid at the inlet would presumably be 



"ideal" pumps are isentropic, but "real" pumps are not. Assuming that 
the efficiency of a "real" pump would be of the order of 50%, we have: 



Let the distance (L-z q ) be designated L', and the slot factor computed with 
L' be S'. Then, 



Actually, the slot factor derived from the previous models is, correctly, 
S'. If S' were multiplied by L, vice L', too large a value would be com- 
puted for s. Presuming L would be known, while L' might not be, the 
previously computed slot factor must be corrected: 



Eqns. (59-61) can be applied to correct the slot factor predicted from 
either flow model. 

Next, the mean temperature increment predicted from the preceding 

models must be re-examined. The original T , which will now be called 

c 

T ', was based on the nonlinear, monotonically decreasing temperature 
profile between z = z q and z = L. By contrast, as discussed previously. 



subcooled by an amount (5h, where 




pump 



(59) 




(60) 



Then, assuming constant (q/A), as before. 




z 



o 



/L ' = <5 h <Pi v i ) TqlWL ' <5h£ 





( 61 ) 
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conditions in the subcooled flow can be approximated as a linear temper- 



ature increase between z=0 and z = z q . (Note that the "hot spot" in the 
windings of a coil cooled by such a system would occur just at z=z q , where 
coolant temperature would roughly equal the coolant saturation temperature, 
at the inlet pressure. ) Since a pump efficiency of 50% would imply rela- 
tively little temperature increase for the coolants considered, even when 
pumped through the maximum pressure increment investigated, the cor- 
rected AT c can be approximated as: 



L'T ' + z (T + T )/2 

o s, pj s, po _ T 

c T s, 1 



AT ~ 



T ' +(5hn.v.S'(T T n )/2 

C ^ r I 1 S,Pj^_+ s,Pe 



- T 



s, 1 



(62) 



1 + ^hpjV.S' 

Finally, it was originally stated that "inlet pressure will be system- 
atically incremented" by assuming p. = p Q + n Ap. Again, the pressures so 
specified correspond to conditions at z = z Q , vice z = 0. With single-phase 
flow to z=z q , Ap could be corrected by deriving an additive term from eqn. 
(30). However, in general z q << L, and we choose to remain with an 
inviscid model. One correction can be applied, though, even to an inviscid 
model: The liquid density decreases between z = 0 and z = z q , since liquid 
density is a significant function of temperature, if not pressure. This 
requires a slight flow acceleration, which results in a pressure decrement 
of the order: 

' 2/i/ ^s r. 

l r o 



(63) 



^inlet-mX^Ps.^-^Ps.p^ 

thus, a "corrected" Ap can be determined by adding eqn. (63) to the Ap 

specified in the incremental analyses. This has been done, though in most 

cases, the correction is insignificant. 

The details of the calculations performed, in applying the no-slip and 

no-mixing models, corrected for subcooling at the inlet, to the analysis of. 
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the two-phase cryogenic cooling problem, are discussed in the next chapter. 
The results of these calculations are presented in Appendix D. 

E. Other coolant flow considerations : 

In closing this chapter, two factors which cannot be discussed quanti- 
tatively deserve qualitative mention: 

The first is the problem of inlet and outlet connections to the coolant 
ducts: although the slot geometry discussed herein probably corresponds 
reasonably well to such geometries as might be employed throughout the 
greater portion of "real" cryogenic saddle coils, the applicability obviously 
terminates at coolant inlet and outlet connections. In fact, no simple 
solution is seen to the problem of compatability with coil geometry at 
these points. Presumably, ducts would not be "spiraled" through several 
axial passes, to emerge at the surface of the coil, since the analysis 
indicates that increased duct length requires increased slot height, and 
therefore decreased coil packing factor. Further, the geometry of the 
"end turns" of the coil is not particularly compatible with conventional 
use of inlet and outlet headers. The problem is a moot one, and regretably 
beyond the scope of this study. 

The second problem has to do with the possibility of a parallel flow 
instability in two-phase cooling systems. It is probable that whatever 
connection scheme were devised would feed several ducts "in parallel". 

Let us assume that steady- state operation has been established, and investi- 
gate what might happen if some local phenomenon perturbs the flow rate 
through a particular duct. 

Fig. 21a depicts two curves of flow rate vs. Ap, such as might be 
characteristic for a particular coolant and a particular geometry at con - 
stant (q/A ). Basically, characteristics of class I are "safe", but those 
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FIGURE 22a 



"TYPICAL*' FLOW CHARACTERISTICS 




FIGURE 22b 

POSSIBLE FLOW INSTABILITY 
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of class II may lead to problems: 



If a high-capacity pump is supplying several ducts in parallel, the 
pump characteristic might imply essentially constant Ap, regardless of 
flow conditions in any one duct. Such a characteristic is shown super- 
imposed on a class II characteristic in Fig. 21b. If steady-state operation 
had been est ablished at point a , and a sudden fluctuation caused a momen- 
tary decrease in mass flow rate in a particular duct, the operating point 
for that duct might "jump" to point b. If the quality at point b were so 
high that flow conditions entered the "liquid deficient" regime, heat trans- 
fer to that duct would then drop abruptly, leading to the creation of a "hot 
spot" extending the axial length of the coil. The problem is obvious; the 
situation can be avoided or ameliorated by careful design; but that is 
also beyond the scope of this study. 
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VII. GENERATION AND USE OF CHARACTERISTIC CURVES 



A. Cur ves of coil geometry : 

Generation : The curves of coil geometry, presented in Appendix A, and 
those of mean coil resistivity, in Appendix B, were generated from a single 
computer program, consisting of a calling sequence (main program), and 
several functional subroutines. A brief resume will be given here of the 
logic used in these routines, and of the numerical ranges investigated dur- 
ing the processing of data. The programs themselves are presented in 
Appendix E. All are written in Fortran, for use with the M. I. T. Compu- 
tation Center version of the Fortran Moniter System, and were run on the 
Computation Center's I. B. M. 7090 digital computer. 

MAIN , as mentioned above, is primarily a calling sequence. As per- 
tains to the coil geometry problem, MAIN reads an input data card, to 
determine the desired channel area (in square meters), desired channel 
flux density (in kilogauss), and effective coil current density (in amps/ 
square meter). The calling sequence is "looped", so that data from one 
card are processed completely, then another card read, etc. As processed 

here, channel area was varied through the specific values of 10, 1, 0. 1, and 
2 

0. 01 m , flux density through 25, 50, 75, and 100 kilogauss, and current 
density varied as necessary to investigate form factors in the approximate 
range 1. 3 = F = 5 (e. g. , 10® to 10^ a/m^ for A ^ = 10, and 5 x 10® to 
10 a/m for A ^ = 0. 01). Customarily, B was varied through all its val- 
ues, for a fixed J g and A then J g incremented, and the process repeated, 
until all coils suitable for a particular channel size had been investigated. 

SIZE performs the actual calculation of coil geometry for a specified 
set of input parameters. Eqns. (1-6) have been programmed in the ration- 
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alized MKS system of units; thus, the specified B is divided by ten to con- 
vert to a consistent set of units. SIZE uses an iterative approach to deter- 
mine the required cylinder radius for a particular coil: eqn. (1) is solved 
for cylinder offset, and a "first guess" is made that 

R ~a + V A ch ’ /2 ' 

This value is inserted in eqns. (2) and (3), and the resulting computed chan- 
nel area compared to the specified channel area. Based on this comparison, 
the estimated R is varied systematically, until the deviation between spec- 
ified and computed A ch is less than 0. 01%, or until the process has been 
repeated 50 times. (The latter option was found to be an unnecessary pre- 
caution, within the range investigated. ) The subroutine then computes the 
width, height, and area of the "optimum" channel for the resultant coil, the 
coil length, total cross section (both lobes), total surface area, and total 
volume (uncorrected for end turns), and finally, the coil form factor. This 
information is stored until MAIN calls the output routine, PRINT. 

Use of coil geometry c urves: The output of the above subroutine is 
presented in Figs. A-l to A- 3, as an ex plicit function of the form factor, 

F, defined in eqn. (6). F has been chosen as the independent variable, for 
two reasons: It provides a convenient, quantitative index to the "shape" 
of a saddle coil, and more importantly, provides an immediate indication 
of design optimization. (This latter point will be discussed in Chapter 
VIII. ) To utilize the information presented in these curves: 

1. Specify desired A^, B, and F. 

2. Enter Fig. A-l with F, and determine required to attain that F. 

3. Enter Fig. A-2 with F. Compute w^, a, and R by multiplying 
plotted values by 10~\J A^ . 

4. Enter Fig. A- 3 with F. Compute A g and A^ by multiplying plotted 
values by 100 A^. Compute V by multiplying plotted value by 1000(A c ^)^^. 
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5. Correct for end turns by multiplying A and V by F. 

s 

As stated in Chapter I, the above assumes that the coil length, L, is 
equal to 10 vA^. However, the information plotted in these curves can be 
corrected for L j lO^A^: 

1. Carry out steps 1-4, as above. 

2. Correct A and V by multiplying by L/ 10 A , . 

s cn 

3. Correct F: F' = 1 + (F - 1)( 10 V A ch 7U. 

4. Correct for end turns by multiplying A and V by F'. 

s 

Note that the optimization procedure discussed in the next chapter per- 
tains to F', vice F, when L is not equal to 10"^A~^. However, all parameters 
plotted as functions of F must be "looked up" as functions of F. Of the geo- 
metric parameters, only A and V need be corrected for the "new" length. 

s 

B. Mean coil resistivity curves: 

Generation : Mean coil resistivity was computed as a function of tem- 
perature, at ten degree increments, from 5° to 95° K, for each coil proc- 
essed by MAIN. To reduce input-output time, the necessary field-free and 
magnetoresistance data for a particular metal were "read in" only once, at 
the start of each run, by having MAIN call READ before entering the com- 
puting loop. 

READ is strictly an input routine, designed to read, in sequence: 

1. Ten specified temperatures, from a single card; 

2. The field-free resistance corresponding to each specified temperature 
(see Tables 1-4), from two data cards; 

3. The magnetoresistance for the particular metal, as tabulated in Table 
5, from ten data cards. 

FIELD determines the field distribution within the upper-right-hand 
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quadrant of the coil, to permit calculation of magnetoresistance. (A modified 
version of FIELD was also used to compute the vector diagrams presented 
in Figs. 12-16. ) The upper-right-hand quadrant was selected, since the 
symmetry of the field distribution indicates that the mean resistivity of any 
quadrant is equal to the mean resistivity of the entire coil. 

FIELD positions a "grid" of 5000 squares over the upper-right-hand 
quadrant of the coil. The overall size of the grid is 2R wide by R high, with 
the lower left corner of the grid positioned at the center of the MHD channel. 
Each component square has sides of dimension R/50. 

FIELD first scans the left-most vertical column of squares, working 
from bottom to top, and determines which squares have centers which fall 
within the coil quadrant. FIELD computes and saves the value of the field 
intensity at the center of these squares, using the equations: 



where the x-y origin is the lower left corner of the grid. After locating 
the uppermost square which falls within the coil quadrant, FIELD shifts to 
the next vertical column to the right, etc. , until locating the first such 
column which lies entirely outside the coil. The field intensity at the center 
of each square which fell within the coil is saved in an array, with the loca- 
tion of the data corresponding to the location of the square in the grid. Loca- 
tions corresponding to squares which fell outside the coil are set to zero. 
FIELD also saves the location of the first column which fell entirely out- 
side the coil. 




( 66 ) 



(64) 



(65) 
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RHO next computes mean coil resistivity, using the field distribution 



generated by FIELD, and the resistance data supplied by READ. RHO 
sequentially examines each specified temperature. For each temperature, 
RHO determines the field-free resistivity, then scans the array of field 
intensities for non -zero entries (examining only those columns which were 
previously scanned by FIELD). RHO divides each non-zero entry by the 
field-free resistivity, and used the result as an index to find the "nearest" 
corresponding entry in the magnetoresistance table. RHO then computes 
"total" resistivity at the center of the corresponding square, and assumes 
this to be representative of the total resistivity of the entire square. Fi- 
nally, RHO computes mean coil resistivity, for the particular temperature, 
from the relation: N 

(67 > 

i= 1 

N being the number of locations found to have non-zero entries. 

In this incremental analysis (as with all other incremental analyses 
conducted herein!), the "fineness" of the grid was selected through suc- 
cessive reduction of the increment used, so that the degree of approximation 
associated with a particular increment size could be estimated. The incre- 
ment specified ensures that the accuracy of the computational technique is 
consistent with the accuracy of the physical model. 

Use of resistivity curves: The data in Tables 1-5 indicate that indium 

coils would not be competitive with coils wound from the other conductor 

materials investigated, except at temperatures below 4° K, v As it seems 

unlikely that cryogenic coils would be designed for operation below 4° K, 

indium was dropped from the list of potential conductor materials, and no 

mean coil resistivity curves are presented for that metal. However, curves 

of mean resistivity vs. temperature for coils wound of sodium, copper, and 
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aluminum, are presented in Figs. B-l to B-16. These curves have been 
normalized to the 1. 7715 micro-ohm-centimeter standard, to permit direct 
comparison with "conventional" copper coils. Further, resistivity curves 
for all three metals are superimposed in each figure, to permit direct com- 
parison between these metals as cryogenic conductors. Figs. B-l to B-4 
are for 25 kg coils, B-5 to B-8 for 50 kg coils, etc. , with each figure pre- 
senting data for a particular form factor. (F and B together determine the 
geometric distribution of magnetoresistance, and therefore define mean 
resistivity, for a particular metal. However, the variation is sufficiently 

slow except for copper below 30° K that interpolation is generally 

unnecessary: to determine mean resistivity, simply enter that figure 
which most nearly corresponds to the contemplated application. 

While these curves were computed from the data presented in Tables 
1-5, approximate adjustment can be made, should it be desired to assume 
some value for any one of the three components of "total" resistivity, 
other than those used here: the plotted "mean" resistivity can be compared 
to the tabulated "field-free" resistivity, to determine the percentage contri - 
bution of each component to total resistance. (Note that the 0° K entry in 
each of Tables 1-4 is the assumed residual resistance ratio for the various 
metals. ) Once the percentage contribution of a component is known, the 
net effect of changing that component can readily be estimated. 

C. Structural curves : 

Generation : The curves of structural volume fraction, presented in 
Appendix C, were generated by a second computer program, consisting of 
a calling sequence (main program), and the subroutines SIZE, CSHELL, 
and TSHELL. 
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MAIN is again "looped"; uses the same input data as the previous pro- 
gram (A B, and J g ), and again calls SIZE to evaluate the geometric 
parameters of each saddle coil so specified. MAIN then systematically 
varies the specified conductor stress limit ( 0 ) through the values 3, 000, 

6, 000, 9, 000, and 12, 000 p. s. i. , and the structural stress limit ( 0^) 
through the values 120, 000, 90, 000, 60, 000 and 30, 000 p. s. i. (All values 
of (JVj are "matched" with a single value of 0^, before the latter is incre- 
mented. ) CSHELL and TSHELL are called for each combination of and 
0^, to compute the compressive and tensile volume fractions, respectively. 

To explain this particular sequencing: The required number of struc- 
tural shells decreases as (5 increases. Thus, if no shells are required 
for a particular 0^, none will be required for higher (f . Accordingly, 
if total computed volume fraction is found to be identically zero for any 
value of O MAIN does not search to higher values, but concludes that 
the problem has been "solved", for that particular coil, and calls for new in- 
put data. 

Similarly, shell thickness increases with decreasing 0^. Thus, if the 
structural volume fraction is unreasonably high, for a given 0^, MAIN does 
not search to lower 0^, but increments O' immediately, to reduce the 
number of shells. "Unreasonably high" is defined, in the program, as a 
volume fraction of 0. 8 or greater. 

SIZE, aside from a minor reduction in the number of "output" parameters 
computed, is the same subroutine as described above. 

CSHELL computes the compressive volume fraction by the procedure 
described in Section D of Chapter III. Since electrical variables and coil 
geometry data are expressed in metric units, while Windenburg's equations 
(and most available strength data) are in English units, considerable con- 
version of units occurs during the computational process. 
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TSHELL computes the tensile volume fraction by the process described 

in Section E of Chapter III. Specifically, TSHELL slices the coil lobe into 

shells of thickness dr = a/ 10, where a is the offset between the center of the 

cylinder and the center of the channel. Per unit length, the volume of each 

such shell is assumed to be 2r9 dr, where r is the inner radius of the shell, 

and 9 is defined by eqn. (20), in terms of that radius. TSHELL integrates 

eqn. (18) numerically, by "piling up 1 such shells until the effective pressure 

exceeds . TSHELL then "shaves" the last shell, reducing its thickness 

so that p = (5 at the outer radius of the stack. A structural support shell 

c 

is inserted at this radius, and its volume determined as indicated in eqn. 
(22). TSHELL continues this process until the outer edge of the coil is 
reached, then divides the accumulated volume by the volume of a unit length 
of coil lobe. 

When the above computation has been completed, TSHELL prints 0 , 
Cf d , the number of compressive and tensile shells, the compressive and 
tensile volume fractions, and the total structural volume fraction. Control 
is then returned to MAIN, which specifies a new combination of O c and 
0^, or begins work on the next coil, as appropriate. 

Use of structural curves : Curves of structural volume fraction as a 
function of form factor are presented in Figs. C-l to C-4, with each figure 
presenting the data for a particular 0^. Interestingly enough, for a given 
0 Q j, structural volume fraction was found to be much more sensitive to F 
and B, than to 0 . Accordingly, curves are presented for each value of 
B investigated, but only for 0 ' - 3, 000 and 12, 000 p. s. i. The structural 
volume fraction for intermediated values of 0 ' can be estimated, to within 
the probable accuracy of the model, by simple linear interpolation between 
these two families of curves. 

D. Coolant curves: 
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Generation: The coolant curves, presented in Appendix D, were gen- 



erated in two stages: First, curve fitting routines were utilized to develope 
the saturation property polynomials, and then, using these polynomials, 
curves of slot factor (S), mean coolant temperature increment ( AT c >, and , 
coolant inlet velocity (v. ) were generated by the techniques discussed in 
Chapter VI. 

The curve fitting routines were executive programs built around an 
M. I. T. Computation Center SHARE program, GLSQ. Briefly, GLSQ deter- 
mines a least-square solution to an overspecified problem (M equations in 
N unknowns, with M > N), or the exact solution to a "properly set" problem 
(M equations in M unknowns), by matrix inversion techniques. 

Both GLSQ, and an example of the type of data fitting routine used in 
conjunction with it, are presented in Appendix E. The data fitting routine 

specifies the functional form to which the data are to be reduced e. g. , 

N- 1 

y = a^ + a^x +---- + a^x and for each data point (x, y), computes 

the various powers of x. (Actually, the functional form need not be con- 
strained to a simple power series in x: it could as well involve terms in 
e , sin x, etc. ) These factors are sent to GLSQ in an array, with each 
row of the array corresponding to one data point. (Obviously, there must 
be at least as many data points as there are coefficients to be evaluated in 
the set {a^} . ) GLSQ solves the matrix for the unknown coefficients {a^ . 

The example program prints these coefficients, the resultant minimum 

mean square deviation, and lists x, 7 observed , y gen erated- ^obs ' ^gen' 

and 100(y o k g -y en )/y en deviation) for each data point. (Input required: 
o s gen gen 

one card, listing the number of data points (M), and M cards, each listing 
the coordinates of a single data point. ) Of the two figures which appear 
repeatedly in the DIMENSION statement, 200 is the maximum number of 
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data points which can be processed, and 7 is the number of terms in one row 
of the matrix. Both of these figures may be adjusted, as desired. 

The coolant flow program consists of a main program, the subroutines 
READS, NOSLIP, NOMIX, INLET, and PRINTS, and the functions SVOLL, 
SVOLV, TEMP, HLIQ, HTVP, AVSVL, and AVSVV: 

READS determines the saturation property coefficients (as listed in 
Table 6) for the coolant to be investigated, the three values of outlet quality 
to be examined, and the limiting velocity to be applied to the flow. 

NOSLIP reads in the outlet pressure to be assumed, the "step size" 
to be utilized in incrementing the inlet pressure, and the number of incre- 
ments to be examined. NOSLIP then solves the no-slip flow model, exactly 
as discussed in the preceding chapter, using a 25-station incremental model. 
Iterative solutions are cycled until the one-cycle change in the control vari- 
able is less than 0. 1%. For a given p. and p Q , flow conditions are solved 
for each specified x q . After this, p. is incremented, and the process re- 
peated until all specified inlet pressures have been investigated. 

NOMIX solves the no-mixing flow model in the same manner. One 
additional point should be noted: Statements 9 and 10 of this subroutine 
generate a corrected estimate for the constant K, which appears in eqns. 
(51-53), in order to obtain convergence between computed and desired out- 
let qualities. Actually, the "loop gain" in this particular iterative solution 
is a strong function of the degree to which specified operating pressures 
approach the critical pressure of the coolant being investigated. The atten- 
uation factors of 10. 0 and 5. 0, which appear in statements 9 and 10, as 
listed in Appendix E, provide suitable convergence for nearly all runs 
covered here. However, if specified pressures lie quite close to the critical 
pressure e. g. , above about 7. 5 atmospheres for liquid hydrogen, or 
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cThe atmosphere and above for liquid helium loop gain is so small that 

the solution converges too slowly if these factors are used. They should 
be replaced by the factors 3. 0 and 1. 5, respectively. (Conversely, if the 
latter values were used for, say, liquid nitrogen in the pressure range of 
1 to 10 atmospheres, the resultant correction would be excessive, and the 
solution would oscillate, or diverge. ) 

INLET corrects the flow solution for subcooling at the inlet, as indicated 
in eqns. (61-63). 

PRINTS can be used to list the output data from all three of the preceding 
analytical subroutines, or, by omission of statements 1-25, can be used to 
obtain only the "corrected" data, generated by INLET. The latter option is 
obviously less time consuming. 

The functions SVOLL and SVOLV utilize the saturation property coef- 
ficients, obtained from READS, to compute, respectively, the specific vol- 
ume of the saturated liquid, and specific volume of the saturated vapor, in 
accord with the functional forms listed in Section B of Chapter VI. Similarly, 
TEMP computes saturation temperature, HLIQ the specific enthalpy of the 
saturated liquid, and HTVP the heat of vaporization, for a specified pres- 
sure. AVSVL and AVSVV respectively compute the mean specific volume 
of liquid and vapor, for use in eqns. (47) and (48), for a specified p. and p Q . 

Use of coolant curves: Curves of S, AT , and v. vs. Ap for each 

'em 

of the various coolants studied are presented in Figs. D-la to D-12c: D-l 
to D-3 pertaining to liquid nitrogen, D-4 to D-6 to neon, D-7 to D-9 to 
parahydrogen, and D-10 to D-l 2 to helium. 

Each figure presents a family of curves, showing variation of a par- 
ticular property with Ap, for various outlet pressures, at a specified 
quality: e. g. , D-la presents data for x q = 30%, D-lb for x q = 50%, and D-lc 

for x = 70%. The "a-b-c" sequence is similarly repeated in all other fig- • 
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ures. Similarly, the D-l set are curves of S vs. Ap, D-2, AT c vs. Ap, 
and D-3, v^ vs. A p. This sequence is also repeated for each coolant. 
Finally, the data generated by the no-slip and no-mixing flow models are 
displayed side-by-side, in each figure, for ease of comparison. 

In general, the plotted curves should not be extrapolated to higher values 
of Ap. Each curve is terminated, either because v , would exceed the 
specified limiting velocity, or because p^ n would exceed the specified limit- 
ing pressure, were Ap increased further. However, the curves could 
certainly be extrapolated to lower values of Ap than those plotted. 

For a specified p Q and Ap, the inlet velocity curves permit estimation 
of the required coolant pump power, as indicated by eqn. (34). Similarly, 
the AT c curves permit estimation of that parameter, for insertion into 
eqn. (58), and the slot factor curves permit estimation of the coolant volume 
fraction: given a preliminary estimate of q (total I R loss), q/A, and L, 

coolant volume = V c = As = ^q/(a/ A)j Js(q/ A)Lj = qSL (68) 

and 

C. V. F. = V c /V , (69) 

where V is as given by eqn. (7). 

Thus, these curves supply the final link in the chain of information 
required for preliminary analysis of cryogenic saddle coils. A suggested 
procedure for conducting such an analysis is discussed in the next chapter. 
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VIII. DESIGN OF CRYOGENIC SADDLE COILS 



A. Ge neral design characteristics : 

Presumably, the first parameters which would be specified in the design 
of an MHD saddle coil would be the desired channel size, and desired field 
intensity: these requirements are set by the intended MHD application, and 
could not conveniently be altered to suit the whim of the coil designer. Since 
we have assumed L = 10 VaC- in line with Sutton's work (28, 29), this implies 
specification of B and A as the first step in the design process. 

Assuming this has been done, the next step in the design of cryogenic 
saddle coils would be the selection, of a particular metal for the coil con- 
ductors, and the selection of a particular coolant (and, by implication, spec- 
ification of the approximate range of possible operating temperatures). "Next", 
because the selection is obviously limited; the choice can be simply and di- 
rectly related to the overall design philosophy, and will have a profound 
influence on subsequent development of the design. 

Figure 23 depicts the total electric and refrigeration power requirement 
(neglecting heat leak and coolant pump power) as a function of coil operating 
temperature, for a "typical" cryogenic saddle coil (B= 50 kg, A ^ = 0. 1 m , 

F = 2), assuming unity packing factor. The approximate range of operating 
temperatures attainable, with each of the liquid coolants examined here, is 
indicated at the bottom of the figure. Although the specified assumptions 
indicate these curves to be a "rough" approximation to the characteristics 
of "real" coils, several conclusions nevertheless appear valid: 

1. Under the same assumptions, a conventional, copper saddle coil 
would dissipate 16. 7 Mw. Comparing this to the curves in Fig. 23, it is 

apparent that liquid nitrogen would not be used to cool "working" saddle coils. 
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2. If liquid nitrogen were used to cool an experimental coil, copper 
conductors would yield lowest power consumption. If coil weight were an 
important factor, aluminum conductors would seem the best compromise 
between weight and power consumption. 

3. In the liquid neon range and below, aluminum appears to be the 
"optimum" conductor material. Optimum operating temperature for aluminum 

saddle coils is about 16° K i. e., about the same temperature recommended 

by Taylor and Post (26) for aluminum solenoids and the optimum is fairly 

sharp: a 10° temperature shift in either direction increases losses by more 
than 50%. 

4. Based on this crude model, operation at liquid neon or liquid helium 
temperatures would appear roughly equivalent, from the viewpoint of power 
consumption. However, a comparison of the slot factors presented in Figs. 
D-4 and D-10 indicates that use of the latter coolant would result in lower 
packing factors, and thus greater losses. Again, liquid helium probably 
would not be used as a coolant for cryogenic saddle coils, except when 
cryogenic conductors are used to "stabilize" superconductors, as discussed 
by Kantrowitz and Stekly (14). 

5. If a two-phase coolant is to be used, liquid hydrogen would seem the 
best choice. 

6. Supercritical helium is the only coolant which could maintain the 
"optimum" temperature for aluminum coils. 

Thus, assuming the motivation for utilizing cryogenic saddle coils is 
that of reduced operating power requirements, one would normally select 
aluminum for the coil conductors, and liquid hydrogen or supercritical helium 
for the coolant, based on these data. 

The next fundamental variable to be considered, after the specification 
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of B, A ch , the conductor metal, and the coolant, is the general "shape" of 

the contemplated coil. While the specified saddle coil model seemingly 

determines this, a glance at Figs. 12-16 will indicate the wide variation of 

form still possible within this family of coils. Of course, "shape" can be 

quantitatively described by the value of the form factor, F, and Fig. 24 

displays relative power consumption as a function of F (under the same 

assumptions that applied to Fig. 2 3). Here, the minimum occurs at F = 

2. 0. Power consumption rises gradually as F is increased above this value: 

2 

loss density (p J 0 ) varies inversely with F, but a comparison of Figs. A-l 
and A- 3 indicates that Je decreases relatively slowly with F, above F = 2, 
while coil volume increases quite rapidly. As a result, total power con- 
sumption and coil size both increa se if F increases above 2. 0. One can 
therefore immediately class form factors above that value as "not of interest". 
Conversely, total power consumption rises if F is decreased below 2, be- 
cause the sharp rise in J more than offsets the decrease in total volume. 

e 

However, though Fig. 24 indicates that a coil with F = 1. 34 would have twice 
the power consumption of one with F = 2, Fig. A-3 indicates that its volume 
would be smaller by more than an order of magnitude. Thus, for some appli- 
cations, form factors less than 2 might be desirable. 

As defined here, the minimum attainable F is theoretically 1. 22: this 
corresponds to a sinusoidally distributed current sheet (infinite current den- 
sity), flowing axially in a cylindrical shell of zero thickness. Realistically, 
Figs. C-l and C-4 indicate that structural requirements alone would con- 
strain "real" coils to form factors above about 1. 3-1. 6, depending upon the 
specified flux density, and allowable stress levels. Thus, w e can define the 
"range of interest" of coil shapes at those corresponding to form factors 



between 1, 3 and 2. 0 , with the lower end of the attainable range corresponding 
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to minimum power consumption. 



B. Simplified design procedure: 

Bearing the above in mind. Fig. 25 presents a flow diagram of the 
recommended design procedure for cryogenic saddle coils. Although the 
motivation for much of the sequencing is obvious, some explanation of var- 
ous points might be useful: 

The design cycle begins with F=2, since both coolant and structural 
volume fractions are lowest at this point, within the range of forms defined 
above. The design process entails an iterative search on these factors, to 
estimate maximum attainable packing factor for a given F. Packing factor, 
of course, is determined from the relationship: 

X = 1 - (S. V. F. + C. V. F. ). (70) 

As F is decreased, the point will be reached where the sum of the coolant 
and structural volume fractions equals or exceeds unity: i. e. , where no 

physical solution exists. However, if a physical solution exists at all 

that is, if it is possible to build any "real" coil, for a given B, A con- 
ductor, and coolant combination the form factor of 2. 0 will correspond 

to a "realizable" coil. Successive reduction of F can then be used to deter- 
mine the smallest "realizable" coil, should that be desired. 

Once B and F have been specified, the structural volume fraction can 
be determined immediately (from Appendix C), since it is not a function of 
resistivity, temperature, coolant pressure, etc. However, to determine 
the required coolant volume fraction, total power dissipation must be known, 
and to determine power dissipation, packing factor must be known. For non- 
unity packing factors, loss density becomes pj = p(J Q / X ) > and total con- 
ductor volume becomes ^FV. Thus, 
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SIMPLIFIED FLOW DIAGRAM FOR SADDLE COIL DESIGN 
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eqn. (58), Chapter 5, and Appendix D 7. Use eqn. (70) 









( 71 ) 



p e (A>= pvFJ e 2 /A =p e (A=D/A 

It is suggested that a "first guess" be made that X = 0. 75 - S. V. F. For 

i 

the smaller F=2 coils, this is reasonably accurate, and for the larger F=2 
coils, it will be found to be conservative. P can then be estimated, the 
corresponding coolant volume fraction computed, and X recomputed from 
eqn. (70). Note that if the "recomputed" X equals or exceeds the "assumed" 
X , the "assumed" X is_physically realizable. However, unless the two 
are equal, the "realizable" coil has not been optimized, in the sense of 
obtaining minimum power consumption for the specified form factor. Thus, 
the design should be recycled until maximum attainable packing factor has 
been determined. (The solution can be made to converge rapidly. In fact, 
after some practice manipulating Appendices C and D, a better "first guess" 
can generally be made than that suggested. ) 

In order to determine coolant volume fraction (and, for that matter, 
mean coil operating temperature) for a two-phase coolant, one must specify 
a particular combination of x q , p Q , and Ap. Presumably, the simplest 
procedure here would be to estimate the highest outlet quality which can 
safely be maintained (since coolant volume fraction decreases as quality 
increases), select a "reasonable" A p (which can be related to coolant pump 
power through eqn. (34) ), and "search" on p Q , starting at one atmosphere 
and working upward. Coolant volume fraction generally (but not always) 
decreases with increased outlet pressure, due to increased vapor density 
in the coolant stream. The question here is whether or not the increased 
packing factor presumably attainable with higher outlet pressures decreases 
power consumption more than the higher coil operating temperature (due to 
higher AT £ ) increases power consumption. (Or, if not, and minimum size 
is being sought, if the higher loss may be tolerable, in exchange for smaller 
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